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FOREWORD 


This  report  presents  a  method  for  calculating  the  inflation 
reference  time  and  opening  shock  forces  of  the  solid-cloth  family 
of  parachutes  when  deployed  at  an  arbitrary  trajectory  angle  to 
the  horizontal . 

The  method  is  extended  to  other  types  of  parachutes,  but  it 
is  limited  in  that  the  inflation  reference  times  are  not 
calculated  as  part  of  the  deployment  process.  Particular 
inflation  times  must  be  provided  as  input  data  to  the  furnished 
computer  program.  The  variation  of  opening  shock  force  versus 
inflation  reference  time  may  be  surveyed  by  providing  several 
values  or  using  actual  field  test  data. 

Examples  are  used  to  demonstrate  the  effects  of  canopy  cloth 
rate  of  airflow,  altitude,  and  trajectory  deployment  angle  for 
constant  velocity  and  constant  dynamic  pressure-altitude  profiles. 

The  author  wishes  to  express  his  appreciation  to  Mr.  Hensel 
Brown  of  the  Strategic  Systems  Department  for  his  valuable 
assistance.  Mr.  Brown  developed  the  computer  programs  and 
conducted  numerous  background  calculations  for  this  report. 
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INTRODUCTION 


Studies  of  parachute  opening  shock  force  in  the  horizontal 
and  vertical  deployment  modes  are  particular  solutions  to  a 
difficult  problem.  Often  stores  are  launched  in  loft  or  dive 
modes.  Other  stores  may  be  launched  horizontally,  but  retarder 
deployment  is  delayed  until  such  time  as  neither  the  closed  mode 
horizontal  analysis  or  the  vertical  approach  is  applicable.  In 
this  report  an  opening  shock  force  technique  is  developed  for  a 
generalized  trajectory  angle  at  deployment.  A  method  for 
calculating  the  inflation  time  of  solid-cloth  parachutes  is 
presented  by  way  of  Mode  1  of  a  two-mode  computer  program.  Mode  2 
of  the  computer  program  is  applicable  to  other  types  of  parachutes 
where  the  inflation  reference  time  is  known  or  can  be  assumed  and 
varied. 
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APPROACH 


The  analysis  is  based  upon  the  application  of  Newton's  second 
law  of  motion  to  particle  trajectories  of  inflating  parachutes 
deployed  at  arbitrary  angles  to  the  horizontal  plane.  Equations 
for  the  horizontal  deployment  of  parachutes  with  a  dynamic  drag 
area  signature  of  the  form  cDs/cDsn  =  ( 1-7")  ( t/t  )^+ 7"  was  developed 
in  Reference  1 .  The  same  drag  area  signature  was  incorporated 
into  the  vertical  toward-the-earth  trajectory  of  Reference  2.  The 
exponent  j  determines  the  type  of  parachute.  A  value  of  j=l  is 
indicative  of  geometrically  porous  ribbon  and  ring  slot  types  of 
canopies.  Solid-cloth  parachutes  have  a  value  of  j=6. 


CpS  =  Instantaneous  parachute  drag  area  during  inflation 

CDS0  =  Steady  state  drag  area  of  the  fully  inflated  parachute 

7*  =  The  ratio  of  the  parachute  drag  area  at  t=0  to  the  steady 

state  parachute  drag  area,  T  =  CpS^/CpS^ 


FIGURE  1.  ARBITRARY  TRAJECTORY  ANGLE  POINT  MASS  FORCE  SYSTEM 
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DEVELOPMENT  OF  THE  VELOCITY-TIME  EQUATION  AND  OPENING  SHOCK  FORCE 
With  reference  to  Figure  1. 


IFX  =  max 


ZFy  =  may 


1  ,  WAV. 

-PV  CDS  Cosy  =  7  — 


1  ,  WAV. 

w-  IpV2CDSSiny  =  -^ 


AV»  = 


-  gpVzCDS  Cos  y 
2W 


At  (1) 


A  .  gpV2CDS  Sin  y  A 
AVv  =  l9  "  2W2 - 'lAt 


(2) 


—  t) 


AVX  =1 


v.,„  ^s.(n-Tl(X)',T)cosA 
vst0  X  2W  / 


(3) 


At  AV  =  l  g  - 


gpV2CDS„  (d  -t)  (^y  +  t)  sin  y' 


At 


2W 


ATt  =  0 


LET  M  = 


2W 


(4) 


LET  At  =  TIME  STEP =1^ 
Vx  =  Vs Cos  y 
Vy  =  Vs  Sin  y 
AT  THE  TIME  t 


IAVX 

Vx=VsCosy 


t 

=  I 

t  =  o 


^-v2(<i-t)  ^y+ t)  c°s^ 


V' 


At  (la) 


MVst0 


2  AVy  =  2 
Vy  =  VsSiny  t  =  o 


g  - 


("  -T’  a) -  T) Sin ' 


MVst0 


At  (2a) 


Change  1 
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Where  the  quantity  2W/(pg  V  t  C  Sn)  has  been  shown  in  Reference  1 
to  be  a  Ballistic  Mass  Ratio  °BMR7  scale  factor. 

The  trajectory  velocity,  V: 


v  =  V  vx2  +  vy2 


(5) 


the  trajectory  angle,  7  : 


Tan  y 


(6) 


At  any  instant  during  the  inflation  process  the  shock  factor  has 
been  shown  in  Reference  1  to  be: 


and  me  msraivcaneous  shock  force  is: 


Fj  —  Xi^s  (8) 

where  Fg  is  the  steady-state  drag  force  of  the  fully  open 
parachute  at  the  line  stretch  velocity,  Vc . 

O 

Fs  =  VipVs2CDS0  ,  . 
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DEVELOPMENT  OF  THE  INFLATION  REFERENCE  TIME  EQUATION  FOR  SOLID- 
CLOTH  PARACHUTES 


The  solid-cloth  parachute  inflation  reference  time  "t  "  in 
the  arbitrary  trajectory  angle  mode  can  be  determined  from°the 
basic  mass  flow  equation  used  to  calculate  the  horizontal 
deployment  "t  ".  Following  the  method  of  Section  IV  of  Appendix  A 
for  solid-cloth  parachutes.  Appendix  B  is  included  as  a  guide  to 
assist  in  the  use  of  Appendix  A  in  the  solution  of  particular 
problems.  Refer  to  page  A-16  for  nomenclature  definitions. 


A  program  for  calculating  tQ  for  solid-cloth  parachutes  and 
the  opening  shock  force  profile  during  the  inflation  of  several 
parachute  types  is  provided  in  Table  1.  Equations  (la)  through 
(11)  are  programmed  together  with  the  opening  shock  equations  (7) 
through  (9)  in  FORTRAN  IV  language.  The  included  examples  were 
calculated  via  the  program  using  a  VAX  780  computer.  The  program 
operates  in  two  modes.  Mode  1,  for  solid-cloth  parachutes, 
calculates  the  deployment  inflation  reference  time  t  for  the 
parachute  system  parameters  and  operational  deploymeRt  data,  and 
also  calculates  the  opening  shock  profile  during  inflation.  A 
typical  data  print  out  is  shown  in  Table  2.  It  is  necessary  to 
estimate  an  initial  value  of  "t  ".  The  program  calculates  the 
canopy  volume  for  the  estimated°time  and  compares  the  to 

the  volume  derived  from  the  canopy  geometry.  If  the  calcufaRed 

dm  =  minfiow  -  "i0utf|0W  where  m  refers  to  mass  of  air  flowing 


P  df  =  pVAm  "  pAsP 

Instantaneous  canopy  mouth  area. 

Am  =  Amo(i) 

Instantaneous  pressurized  canopy  cloth  area. 


As  =  As 


(if 


Instantaneous  cloth  rate  of  airflow. 

c_ 


P  =  k 


'P  P\  y2n 


FOR  t  =  0 


FOR  t  =  0 


C.F.S./FT2 


Where  the  instantaneous  velocity  is  determined  by  equation  (5)  for 

r  =  0. 

The  calculated  canopy  volume  iT  „  .  .  is  determined  from  equation 

(10).  ca±c 


ocalc 


V0  t0 

=  I  AY  =  I 
Y=  o  t  =  o 


VAmo  (- 

f)’— (; 

6 

'-)  k< 

;vH 

\ 

o  / 

At 


(10) 
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TABLE  1.  PARACHUTE  PROGRAM 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 


***#**##*#****#**#####*  tt****.****.*#.*.*.*.*.*#.*.*.*.*.*.*#.*#.*.*.*^.*#.*#^^^^^^^^^^^^ 

THIS  PROGRAM  CALCULATES  INSTANTANEOUS  VELOCITY,  DISTANCE  OF  FALL,  ft 
DRAG  AREA,  AND  OPENING  SHOCK  FORCE  OF  A  PARACHUTE  DEPLOYED  IN  ft 

ANY  DIRECTION.  » 

THE  PROGRAM  OPERATES  IN  TWO  MODES. 


MODE  1  -  CALCULATES  THE  INFLATION  TIME  AND  PERFORMANCE  PROFILES  FOR 
SOLID-CLOTH  PARACHUTES  (TO  INPUT  AS  INITIAL  ESTIMATE) 

( I DPT= 1 ) 


MODE  2  -  CALCULATES  THE  PERFORMANCE  PROFILES  FOR  VARIOUS  TYPES  OF 
PARACHUTES  (J).  INFLATION  TIME  IS  REQUIRED  AS  INPUT 
( I 0PT=2 ) 


INPUT  NEEDED  FOR  BOTH  MODES: 


ALT 

RHO 

VS 

CDSO 

TO 

w 

J 

TAU 
ALPHA ( 0 ) 


ALTITUDE.  (FT)  NOT  USED  IN  ANY  CALCULATIONS. 

AIR  DENSITY  AT  GIVEN  ALTITUDE.  ( SLUGS/FT3 ) 

VELOCITY  AT  SUSPENSION  LINE  STRETCH  (FT/SEC) 

DESIGN  DRAG  AREA.  ( FT2 ) 

I OP T= 1  INITIAL  GUESS  FOR  INFLATION  REF  TIME.  (SEC) 
I  OP  T=2  ACTUAL  INFLATION  REFERENCE  TIME  (SEC) 
WEIGHT  (LBS) 

=6  FOR  FLAT  CIRCULAR  PARACHUTE. 

=1  FOR  RIBBON  TYPE  OF  PARACHUTE. 

INITIAL  DRAG-AREA  RATIO. 

INITIAL  DEPLOYMENT  ANGLE  (DEG)  POSITIVE  DOWNWARD 


C 

C 

C 

C 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


ADDITIONAL  INPUT  NEEDED  FOR  IOPT  =  1  ONLY 

AMO  -  STEADY-STATE  MOUTH  AREA  (FT2) 

ASO  -  CANOPY  DESIGN  SURFACE  AREA.  ( FT2 ) 
K  -  CLOTH  PERMEABILITY  CONSTANT 

CP  -  PRESSURE  COEFFICIENT 

N  CLOTH  PERMEABILITY  EXPONENT 

VO  -  GEOMETRIC  VOLUME.  ( FT3 ) 


REAL *4  N 
T  ODEN” 1 OOOOO 
5  PRINT  *, 'INPUT  IOPT ' 

READ ( 5; */ END- 1 00 )  IOPT 

PRINT  *,  '  INPUT  ALT,  RHO,  VS,  CDSO,  TO,  W,  J,  TAU,  ALPHA(O)  ' 

READ  (  5,  *  ) ,  ALT,  RHO,  VS,  CDSO,  TO,  W,  J,  TAU,  ALPHAO 

IF(  IOPT.  EQ.  2)  GO  TO  3 

PRINT  ft,  ' INPUT  AMO,  ASO, K, CP,  N,  VO' 

READ(5,  ft)  AMO,  ASO,  XK,  CP,  N,  VO 
3  DT=TO/TODEN 
G-32  2 
S=0 
X  =  T  AU 
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TABLE  1.  PARACHUTE  PROGRAM  (CONT.) 


CDS=TAU*CDSO 

FS=.  5*RH0*VS**2*CDS0 

F=T  AU-»-FS 

V0L=0 

I P ASS=0 

DT  =10.  *DT 

IF< I OPT.  EQ  2)  GO  TO  8 

7  V^VS 

ALP HA= ALPHAO 
V0L=0 
1  “0 

XM= ( 2.  *W)/(RHO*G*CDSO*VS*TO) 

6  TFAC  =  < 1.  -TAU) *< T/TO )**J  +  TAU 

DVX=-V**2*C0SD( ALPHA )*TFAC/  (  XM*VS*TO  )  *DT 
DVY=  (  G-V**2/  (  XM*VS*TO  )  *TFAC*SIND  ( ALPHA  )  )  *DT 
VX=V*COSD< ALPHA) 

VY=V*SIND< ALPHA) 

VX=VX+DVX 
VY= VY+DVY 

V— SORT  <  VX  -*-«-2+VY*-*-2  ) 

DVOL=  (  V-w-AMO*  ( T/TO )  J-ASO-#-  ( T/TO  )  J-h-XK#  ( CP*RH0/2.  )*#N 
1  *V**  (2.  -a-N  )  )  *DT 
VOL=VOL+DVOL 

IF(VX.NE.  0  )  ALPHA=ATAND ( VY/VX ) 

T=T +DT 

IF(T.  GE.  TO)  GO  TO  4 
GO  TO  6 

4  IF  ( VOL.  GT.  (V0-10.  )  .AND  VOL.  LT.  (V0+10.  ))  GO  TO  8 
TO=TO*< VO/VOL) 

IF ( IP  ASS  GT.  50)  STOP  'VOLUME  ITERATION  FAILED  AFTER  50  PASSES' 

IP ASS= IP ASS+ 1 
GO  TO  7 

8  CONTINUE 

XM=  (  2  *W)  /  (RHO*G-k-CDSO*VS*TO> 

T=0. 

v=vs 

ALPHA=ALPHAO 

R=0. 

S^O. 

DS=0. 

DPR INT=T0/50. 

DT =  DT / 1 0 
WR I TE ( 6,  10) 

10  FORMAT ( 1H1 'DEPLOYMENT  CONDI T I ONS ' ,  1 7X,  'SYSTEM  PARAMETERS') 

WRITE (6, 1 1 )  ALT, CDSO 

11  FORMAT  (  1H0,  3X,  'ALTITUDE  =',F7  0,  '  FT',15X,  'STEADY  STATE  DRAG  ARFA 
1'  CDSO  =  ',F8  2,  '  FT2 '  ) 

WR I TE ( 6,  12)  RHQ,  W 

12  FORMAT (  1H  ,  4X,  'DENSITY  =  ',F10.  7,  '  SLUGS/FT3  ' ,  24X,  'WEIGHT  W  =  ', 

1  F6. 1 , '  LBS ' ) 

WR I TE ( 6,  13)  VS,  TO 

13  FORMAT  (  1 H  ,  'VELOCITY  VS  =  ',F6.  1,  '  FT/SEC ',22X,  'INFLATION  TIME  TO 
1  ',  F7  3,  '  SEC  '  ) 

WR  I  TE  (  6,  14)  J,  FS 

14  FORMAT ( 1H  ,  10X,  'J  =',I2, 24X,  'STEADY  STATE  DRAG  FROCE  FS  = ' ,  F9  1, 

1  '  LBS'  ) 

WR I TE ( 6, 15)  TAU, XM 

15  FORMAT  (  1H  ,  BX,  'TAU  =  ',  F6.  3,  45X,  'M  =',  F8.  5) 

WR I TE ( 6, 16)  ALPHAO 

16  FORMAT  (  1 H  ,  5X,  'ALPHAO  =  ',F6  2,  '  DEG') 

IF(  IOPT.  EQ.  1  )  WRITE (6,  19 )  N 
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TABLE  1.  PARACHUTE  PROGRAM  (CONT.) 

19  FORMA  7  <  1  bH  ,  64X,  'N  =',F6.  3) 

I F ( I OPT  EQ  1)  WRITE<6,  17)  AMO.  ASO,  XK;  VO/  CP 

17  FORMAT  (  62X,  '  AMO  =  S  F8.  2,  '  FT2'/ 

1  62X,  ' ASO  =/,F8.  2,  '  FT2  V 

2  64X,  'K  =  S  F6.  2/ 

3  58X,  'VO  GEOM  =',F9  2.  '  FT3  '  / 

4  61X,  'C.  P  =  '  F6.  2) 

WRITE  <6,  19) 

18  FORMAT ( 1HO ) 

WRITE (6, 20) 

20  FORMAT  ( 2X,  'TIMES  3X>  '  VELOC  I  TV  S  4  X ,  'ALPHAS  4X,  'RANGES  7X»  'SS 

1  6X ,  'DISTANCES  6X,  'SHOCKS  9X,  'F'/ 

2  3X ,  SEC  S  4  X  >  'FT/SECS6X,  'DEGS7X,  'FTSBX,  'FTS8X.  'FTS8X, 

3  'FACTOR  S  8X ;  'LBS  '  ) 

LCQUNT  =11 

F--0. 

I B YP ASS=0 

I F  (  I  OPT.  EQ  1  )  LCOUNT  =16 

C 

C  IBYPASS=0  FORCE  HAS  NOT  REACHED  A  MAXIMUM,  CONTINUE  CALCULATION 
C  I B YPASS=1  FORCE  HAS  REACHED  A  MAXIMUM.  STOP  CALCULATION  IF  TITO 

35  IF(T/TO.  GT.  1  001  AND.  IBYPASS.  NE  0)  GO  TO  5 
IF ( LCOUNT.  LT.  54)  GO  TO  33 
WRITE <6, 37) 

37  FORMAT ( 1 HI ) 

WR  I TE  (  6,  20) 

LC0UNT=3 

38  LCOUNT=L.COUNT+ 1 

WR  I  TE  (  6.  30)  T,  V.  ALPHA,  R.  S,  DS,  X,  F 
30  FOR  MAT  (  IX.  F5.  3.  F10,  2,  F10.  2,  F9  1 , 2F 1  1  2.  1  PE  1  3  3.0PF11  1) 

TOLD-T +DPR I NT 

I F ( DT  GT.  DPR  I NT )  TOLD=T+DT 
40  TPREV-T 
7 -T +  DT 

IF(T.  GT.  TOLD)  T =TOLD 

TOVTO=< 1.  ~TAU)*<T/T0)**J  +  TAU 

CDS=CDS0*T0VT0 

DVX  =  ~V**2*C0SD  (  ALPHA  )  *TOVTO/  (  XM-*VS*T0  )  *DT 
DVY= ( G-V**2*T0VT0*SIND ( ALPHA ) / ( XM*VS*TO ) ) *DT 
VX=V*COSD< ALPHA) 

VY=V*S I ND ( ALPHA ) 

VX=VX+DVX 

VY=VY+DVY 

V=SGRT  (  VX*-h-2+VY-»--h-2  ) 

IF ( VX  NE.  0.  )  ALPHA=ATAND( VY/VX ) 

R=R+VX*DT 

S=S+VY*DT 

DS=DS+V*D7 

X=  ( V/VS)-*-ff-2*CDS/CDS0 

FOLD=F 

F~X-*FS 

IF( IB YPASS.  EQ  1)  GO  TO  99 
IF(T-TPREV.  LT.  DT)  GO  TO  99 
IF ( F  GE. FOLD)  GO  TO  99 
IB YP ASS= 1 

WR  I  TE  ( 6.  98)  T,  V,  ALPHA,  R,  S,  DS,  X,  F 

98  FOR  MAT  (  IX,  F5.  3,  F10  2,F10.  2,  F9.  1,  2F11.  2,  1PE13  3,  0PF11  1 

1  7  MAX  FORCE' )  i# 

LC0UNT=LC0UNT+1 

99  IF ( T. LT. TOLD)  GO  TO  40 
GO  TO  35 

100  STOP 
END 
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TABLE  2.  GENERAL  TRAJECTORY  TECHNIQUE  MODE  1  COMPUTER  SOLUTION  OF  EXAMPLE  1 : 
TOWARD-THE-EARTH-DELIVERY  CONDITION,  7  -  90° 

DEPLOYMENT  CONDITIONS  SYSTEM  PARAMETERS 


ALTITUDE  » 
DENSITY  » 
VELOCITY  VS  * 

J  - 

TAU  - 
ALPHAO  « 


3000.  FT 

0.  0021753  SLUGS/FT3 
325.  3  FT/SEC 
6 

0.  000 

90.  00  DEO 


STEADY- STATE  DRAG  AREA  CD50  - 
WEIGHT  U  * 
INFLATION  TIME  TO  - 
STEADY-STATE  DRAG  FROCE  FS  * 

M  » 
N  » 

AMO  * 
ASO  * 

K  » 

VO  QEOM  - 
C.  P.  * 


725.  5Q  FT2 
200.  0  LBS 
0.  751  SEC  . 
83495.  1  LBS 
0.  03221 
0.  500 

399.  53  FT2 
962.  10  FT2 
1.  46 

4690.  83  FT 3 
1.  70 


time 

VELOCITY 

ALPHA 

SEC 

FT/SEC 

DEG 

0.  000 

325.  25 

90.  00 

oi5 

325.  71 

90,  00 

■:*.  020 

326.  17 

90.  00 

:  045 

326.  63 

90.  00 

:.  060 

327.  oa 

90.  00 

0  075 

327.  54 

90.  00 

0*0 

328.  00 

90.  00 

•:  :05 

328. 46 

90.  00 

:  ;2o 

328.  92 

90.  .00 

0  135 

329  38 

90.  00 

0  130 

329  84 

90.  00 

0  163 

330.  30 

90.  00 

0.  ISO 

330.  75 

90.  00 

0  193 

331.  19 

90.  00 

-■  2 1 0 

331.  59 

90.  00 

0  223 

331.  95 

90.  00 

0.  240 

332.  24 

90.  00 

0  235 

332.  45 

90.  00 

-V  270 

332.  55 

90.  00 

0  2S6 

332.  49 

90.  00 

0  301 

332.  23 

90.  00 

0.  316 

331.  71 

90.  00 

0  331 

330.  86 

90.  00 

*:  246 

329.  61 

90.  00 

36i 

327.  87 

90.  00 

0*  3  7  a 

323.  52 

90.  00 

J  391 

322.  47 

90.  00 

:  406 

318.  60 

90.  00 

0 

313.  79 

90.  00 

:•  436 

307.  94 

90.  00 

:  d.  5  i 

300.  94 

90.  00 

J  4*0 

292  75 

90.  00 

0  43 1 

283.  31 

90.  00 

0.  4*6 

272.  66 

90.  00 

0.  511 

260.  94 

90.  00 

0.  526 

247.  98 

90.  00 

0  34; 

234.  24 

90.  00 

•:  53o 

219.  82 

90.  00 

:  57i 

204  98 

90.  00 

0  591 

194.  93 

90.  00 

5.  58o 

189  95 

90.  00 

0.  -60 1 

174.  99 

90.  00 

0  616 

160.  33 

90.  00 

0  621 

146.  17 

90.  00 

0  646 

132.  68 

90.  00 

0.  s61 

119.  97 

90.  00 

0.  676 

108.  13 

90.  00 

0.  691 

97.  21 

90.  00 

0.  706 

87.  21 

90.  00 

0.  721 

78.  13 

90.  00 

0  736 

69.  92 

90.  00 

0  731 

62.  54 

90.  00 

RANGE 

s 

DISTANCE 

FT 

FT 

FT 

0.  0 

0.  00 

0.  00 

0.  0 

4.  89 

4.  89 

0.  0 

9.  79 

9.  79 

0.  0 

14.  70 

14.  70 

0.  0 

19  61 

19.  61 

0.  0 

24.  53 

24.  53 

0.  0 

29.  46 

29.  46 

0.  0 

34  39 

34.  39 

0,  0 

39.  33 

39.  33 

0.  0 

44.  28 

44.  28 

0.  0 

49.  23 

49.  23 

0.  0 

54.  19 

54.  19 

0.  0 

59.  .16 

59.  16 

0.  0 

64.  13 

64.  13 

0.  0 

69.  12 

69.  12 

0.  0 

74.  10 

74.  10 

0.  0 

79.  09 

79.  09 

0  0 

84.  09 

84.  09 

0.  0 

89.  08 

89.  08 

0  0 

94.  08 

94.  08 

0.  0 

99  07 

99.  07 

0.  0 

104.  06 

104.  06 

0.  0 

109.  04 

109.  04 

0.  0 

114.  01 

114.  01 

0.  0 

118.  95 

118.  95 

0.  0 

123.  86 

123.  86 

0  0 

128.  73 

128  73 

0.  0 

133.  55 

133.  55 

0.  0 

138.  30 

13S.  30 

0  0 

142.  98 

142.  98 

0.  0 

147.  55 

147.  55 

0.  0 

152.  02 

152.  02 

0.  0 

156.  35 

156.  35 

0.  0 

160.  53 

160.  53 

0.  0 

164.  54 

164.  54 

0.  0 

168.  36 

168.  36 

0.  0 

171.  99 

171.  99 

0.  0 

175.  40 

175.  40 

0.  0 

178. 59 

178.  59 

0.  0 

180.  60 

180.  60 

0.  0 

181.  56 

181.  56 

0.  0 

184,  30 

184.  30 

0.  0 

186.  82 

186.  82 

0.  0 

189.  12 

189.  12 

0.  0 

191.  22 

191.  22 

0.  0 

193.  12 

193.  12 

0.  0 

194.  83 

194.  83 

0.  0 

196.  37 

196. 37 

0.  0 

197.  76 

197.  76 

0.  0 

199.  00 

199.  00 

0.  0 

200.  11 

200.  11 

0.  0 

201.  10 

201.  10 

SHOCK 

F 

FACTOR 

LBS 

0. OOOE+OO 

0.  0 

6. 418E-1 1 

0.  0 

4.  1 19E-09 

0.  0 

4.  705E-03 

0.  0 

2.  65 IE— 07 

0.  0 

1.  014E-06 

0.  1 

3.  037E-06 

0.  3 

7. 679E-06 

0.  6 

1. 716E-05 

1.  4 

3.  488E-05 

2.  9 

to.  5S2E-05 

5.  5 

1.  169E-04 

9  8 

1. 976E-04 

16.  5 

3.  203E-04 

26.  7 

5.  009E-04 

41.  3 

7.  593E-04 

63.  4 

1. 120E-03 

93.  5 

1.  ol4E-03 

134.  7 

2.  276E-03 

190  0 

3. 146E-03 

262.  7 

4.  274E-03 

356  8 

5.  709E-03 

476.  6 

7.  509E-03 

626.  9 

9  730E-03 

812.'  3 

1. 243E-02 

1037.  6 

1.  565E-02 

1306.  7 

1.  943E-02 

1622.  5 

2.  379E-02 

1986.  2 

2.  87 IE— 02 

2396  5 

3  412E-02 

2848.  8 

3. 994E-02 

3334  to 

4. 60 IE— 02 

3841  5 

5.  214E-02 

4353  0 

5.  808E-02 

4849. 2 

6/359E-02 

5308. o 

6.  839E-02 

5709.  4 

7.  226E-02 

6032.  2 

7.  501E-02 

6262. 0 

7.  to53E-02 

6389  5 

7.  6S6E-02 

6416  5  H AX  FORCE 

7  toSlE-02 

6412.  6 

7  588E-02 

6335.  1 

7.  387E-02 

6167.  3 

7.  G95E-02 

5923. 4 

6.  732E-02 

5620. 2 

to. 31SE-02 

5274.  9 

5.  S74E— 02 

4904.  0 

5.  416E-02 

4522.  0 

4.  960E-02 

4141.  0 

4.  516E-02 

3770.  5 

4.  093E-02 

3417.  4 

3. 697E-02 

3086.  3 
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volume  is  not  within  specified  limits,  the  program  adjusts  "t  "  by 
equation  (11)  and  reiterates  the  program  until  the  calculated0 
volume  is  within  the  specified  limits.  The  geometric  volume  of 
the  parachute  is  determined  using  the  methods  of  Section  VIII  of 
Appendix  A. 

V0  geometric 


Mode  2  of  the  program  calculates  opening  shock  profiles  for 
input  values  of  t  .  Mode  2  analysis  of  other  types  of  parachutes 
is  possible  by  the  selection  of  the  proper  values  of  "j"  (1/2,  1, 
2,  3,  4,  5,  or  6)  and  "T  ".  In  most  applications  T  may  be 
assumed  to  be  zero. 

When  a  given  system  is  deployed  at  the  same  altitude  and 
velocity,  but  at  various  trajectory  angles,  the  inflation 
reference  time  is  a  maximum  for  the  vertical  away-f rom-the-earth 
trajectory  and  a  minimum  for  the  vertical  toward-the-earth 
trajectory.  The  inflation  reference  times  and  opening  shock 
forces  vary  accordingly  with  trajectory  deployment  angle  between 
these  limits.  As  the  altitude  increases  the  variation  of  opening 
shock  force  with  trajectory  angle  is  minimized. 

Equation  (10)  is  the  identical  expression  as  used  in 
Reference  2  for  the  vertical  toward-the-earth  mode  of  delivery. 

The  velocity  term  of  equation  (10)  is  the  only  element  not 
associated  with  the  parachute  geometry  or  canopy  cloth  airflow. 
x he  inflation  time  is  determined  by  the  trajectory  velocity 
variation,  and  the  velocity  variation  depends  on  the  trajectory 
deployment  angle. 

This  method  is  best  suited  to  the  solid-cloth  family  of 
parachutes.  It  does,  however,  give  a  choice  of  methods  for 
horizontal  deployments.  Reference  1  assumes  that  the  trajectory 
remains  horizontal  throughout  the  inflation  process  and  provides  a 
convenient  closed  form  method  of  calculating  opening  shock  force 
which  does  not  require  the  use  of  a  computer.  This  method 
requires  a  computer,  but  does  include  the  two-dimensional 
trajectory  space  variation. 


EXAMPLES  OF  APPLICATION 

As  a  demonstration  of  the  method,  a  solid-cloth  parachute  of 
35  feet  D  diameter  and  30  gores  is  to  be  used  to  retard  a  200- 
pound  weight.  Determine  the  effects  of  trajectory  deployment 
angle,  altitude,  and  cloth  rate  of  airflow  on  the  inflation  time 
and  opening  shock  force  when  deployed  at  constant  velocity  and 
constant  dynamic  pressure. 
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1 .  Parachute  geometry 

a.  Parachute  surface  area: 


S0  =  l  <35)2 


S0  =  962.1  ft2 

b.  Parachute  drag  area:  Preliminary  design  CD=0.75 


C0So  =  0.75  x  962.1 


CDS0  =  721.6  ft2 


From  Table  2,  page  A-15  of  Appendix  A,  the  following  inflated 
shape  data  were  obtained  for  a  30-gore  solid-cloth  flat  circular 
parachute.  See  Figure  24  of  Appendix  A  for  shape  nomenclature. 

—  =  0.668;  ^  =  0.827;  ^  =  0.6214;  =  0.7806 

D0  a  a  a 

c.  Steady-state  inflated  canopy  radius: 

.668  D0 


0.668  x  35 

a  =  - 

2 


a  =  1 1.69  ft. 


d.  Steady-state  mouth  area: 

-(°-8^7-8°66214)2] 


AMo  =  399.53  ft2 
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e.  Steady-state  canopy  volume  of  air  to  be  collected, 
JZQ  geometric. 


V0 

v0 


2  -3 
3na 


[M] 


|  7t(1 1 ,69)3 


6214  +  0.7807 


Y0  =  4690.83  ft3 


The  program  of  Table  1  was  initiated  to  determine  the 
inflation  reference  time  and  opening  shock  force  profiles  for  the 
example  1  parachute  system  deployed  at  various  trajectory  angles 
from  vertigal  away-f rom-the-earth  (-90°)  to  vertical  toward-the- 
earth  (+90  )  at  increments  of  15  degrees.  The  analysis  is 
conducted  at  altitudes  of  sea  level,  3,000,  and  20,000  feet.  At 
each  altitude  the  effects  of  a  constant  velocity  of  325.25  fps  and 
a  constant  dynamic  pressure  of  125.7  psf  are  examined.  Cloth 
airflow  variations  are  illustrated  by  assigning  the  cloth  airflow 
exponent  "n"  (equation  10)  values  of  0.5  and  0.632.  The 
assumption  of  n=0.5  allows  a  closed  form  solution  for  t  in 
Reference  1.  The  n=0.632  is  a  measured  value  for  the  c?oth  of  the 
illustrative  parachute.  The  results  of  example  1  are  presented  in 
Figures  2  through  5,  and  Tables  2  and  3. 

For  any  given  altitude  at  constant  velocity  the  maximum 
inflation  reference  time.  Figure  2,  minimum  opening  shock  force. 
Figure  3,  and  minimum  time  of  occurrence  in  the  inflation  process. 
Table  3,  occur  in  the  vertical  trajectory  away  from  the  earth. 
Conversely,  the  minimum  inflation  reference  time,  maximum  opening 
shock  force  and  maximum  time  of  occurrence  in  the  inflation 
process  occur  in  the  vertical  trajectory  toward-the-earth.  When 
the  value  of  the  cloth  airflow  parameter, n,  is  varied  from  0.5  to 
0.632,  the  inflation  reference  times  are  extended  with  a 
subsequent  reduction  of  the  maximum  opening  shock  forces  and  time 
of  occurrence  of  the  maximum  force  at  all  trajectory  angles. 

As  the  altitude  increases,  the  inflation  reference  time 
decreases.  Figure  2;  the  maximum  shock  force  increases.  Figure  3; 
and  the  maximum  force  occurs  later  in  the  inflation  process 
Table  3. 

When  a  constant  dynamic  pressure  analysis  profile  is 
substituted  for  the  constant  velocity  profile  the  effects  of 
trajectory  angle,  n,  and  deployment  altitude  show  the  same  trends 
but  different  magnitudes  than  the  constant  velocity  approach.  See 
Figures  4  and  5  and  Table  3.  At  sea  level  the  constant  velocity 
and  constant  dynamic  pressure  results  are  identical.  At  all  other 
altitudes  the  parachute  opens  faster  with  an  increase  in  opening 
shock  force. 

The  variation  of  opening  shock  force  as  a  function  of 
trajectory  launch  angle  tends  to  equalize  as  the  altitude 
increases.  Table  4  illustrates  the  force  variation  for  the 
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CONSTANT  VELOCITY  PROFILE. 
POLAR  PLOT  OF  THE  RESULTS 
OF  EXAMPLE  1 


VELOCITY  =  325.25  fps 

-GO” 


-  30,;I 


45': 


LOFT 


INFLATION  REFERENCE  TIME  t 

o m - 1 — 

1.0  1.5 


0  5 


SECONDS 


DIVE 


4  5 


75* 


ALTITUDE 

ft. 

SYMBOL 

n  =  0.5 

n  =  0.632 

0 

o 

• 

3,000 

□ 

■ 

20,000 

A 

A 

LU 

n  0.500 

n  =  0  662 

ALTITUD 

Y 

(deg) 

*0 

(sec) 

*0 

(sec) 

-90 

1.054 

1.416 

-75 

1.047 

1.401 

-60 

1.026 

1.362 

-45 

0.998 

1.31  1 

-30 

0.970 

1.259 

_ i 

LU 

> 

-15 

0.940 

1.213 

LU 

0 

0.913 

1.171 

< 

LU 

15 

0.894 

1.138 

CO 

30 

0.874 

1.108 

45 

0.860 

1.088 

60 

0.849 

1.074 

75 

0.841 

1 .065 

90 

0.839 

1.062 

-90 

0.904 

1.171 

-75 

0.900 

1.163 

-60 

0.894 

1.142 

-45 

0.876 

1.111 

-30 

0.854 

1.079 

4-* 

4- 

-15 

0.833 

1.046 

o 

o 

0 

0.813 

1.018 

o 

CO 

15 

0.794 

0.992 

30 

0.779 

0.972 

45 

0.767 

0.955 

60 

0.758 

0.944 

75 

0.753 

0.938 

90 

0.751 

0.936 

-90 

0.513 

0.590 

-75 

0.512 

0.589 

-60 

0.509 

0.585 

-45 

0.504 

0.579 

-30 

0.499 

0.572 

+-» 

4— 

o 

-15 

0.493 

0.565 

o 

o 

0 

0.486 

0.559 

o 

CN 

15 

0.481 

0.551 

30 

0.476 

0.545 

45 

0.471 

0.539 

60 

0.469 

0.537 

75 

0.467 

0.534 

90 

0.466 

0.534 

FIGURE  2.  EFFECTS  OF  TRAJECTORY  LAUNCH  ANGLE,  ALTITUDE,  AND  CLOTH  RATE  OF  AIRFLOW 
ON  THE  PARACHUTE  INFLATION  REFERENCE  TIME  AT  CONSTANT  VELOCITY, 

EXAMPLE  1 
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CONSTANT  VELOCITY  PROFILE. 
POLAR  PLOT  OF  THE  RESULTS 
OF  EXAMPLE  1. 


0° 


90° 


75° 


-30° 


Y  =  -90 


75c 


VELOCITY  =  325.25  fps 

—  60  0 
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FIGURE  3.  EFFECTS  OF  TRAJECTORY  LAUNCH  ANGLE,  ALTITUDE,  AND  CLOTH  RATE  OF  AIRFLOW 
ON  THE  MAXIMUM  OPENING  SHOCK  FORCE  AT  CONSTANT  VELOCITY,  EXAMPLE  1 
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CONSTANT  DYNAMIC  PRESSURE 
PROFILE.  POLAR  PLOT  OF 
THE  RESULTS  OF  EXAMPLE  1 . 
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FIGURE  4.  EFFECTS  OF  TRAJECTORY  LAUNCH  ANGLE,  ALTITUDE,  AND  CLOTH  RATE  OF  AIRFLOW 
ON  THE  PARACHUTE  INFLATION  REFERENCE  TIME  AT  CONSTANT  DYNAMIC 
PRESSURE,  EXAMPLE  1 


16 


NSWC  TR  88-6 


ALTITUDE  | 

n  =  0.500 
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FIGURE  5.  EFFECTS  OF  TRAJECTORY  LAUNCH  ANGLE,  ALTITUDE,  AND  CLOTH  RATE  OF  AIRFLOW 
ON  THE  MAXIMUM  OPENING  SHOCK  FORCE  AT  CONSTANT  DYNAMIC  PRESSURE, 
EXAMPLE  1 


17 


MSWC  TR  88-6 


TABLE  3.  EFFECTS  OF  TRAJECTORY  LAUNCH  ANGLE,  ALTITUDE,  AND  CLOTH  RATE  OF 
AIRFLOW  ON  THE  TIME  OF  OCCURRENCE  OF  THE  MAXIMUM  SHOCK  FORCE 
DURING  THE  INFLATION  PROCESS 


CONSTANT 

VELOCITY 

CONSTANT 

DYNAMIC 

PRESSURE 

ALTITUDE 

7 

(DEG) 

n  =0.500 

n  =0.632 

n  =0.500 

n  =0.632 

t/t0 

^max 

t/ 10 

@^max 

t/t0 

^max 

t/t0 

®^max 

90 

.7296 

.6995 

0.7296 

0.6989 

-75 

.7307 

.7002 

0.7309 
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0.8880 
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-75 

.8848 

.8676 

0.8877 

0.8661 
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.8861 

.8684 

0.8874 
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TABLE  4.  EFFECTS  OF  ALTITUDE  AND  CANOPY  CLOTH  RATE  OF  AIR  FLOW  ON 
THE  RATIO  OF  THE  MAXIMUM  OPENING  SHOCK  FORCE  OF  THE  VERTI 
CAL  TOWARD-THE-EARTH,  AND  VERTICAL  AWAY-FROM-THE-EARTH 
TRAJECTORIES  AT  CONSTANT  VELOCITY  AND  DYNAMIC  PRESSURE 


R  = 

Fmax@y  =  90° 

Fmax@V  =  -90° 

ALTITUDE 

CONSTANT  VELOCITY 

CONSTANT  DYNAMIC  PRESSURE 

(ft.) 

V  =  32 

?5.25  fps 

q  =  125.7  psf 

n  =  0.5 

n  =  0.632 

n  =  0.5 

n  -0.632 

0 

1.48 

1.64 

1.48 

1.64 

3K 

1.40 

1.51 

1.37 

1.46 

20K 

1.23 

1.25 

1.11 

1.13 

50K 

1.16 

1.17 

1.021 

1.022 

100K 

1.155 

1.157 

1.0012 

1.0016 

toward-the  earth  and  the  away— f rom- the -earth  trajectories  as 
affected  by  altitude  and  the  rate  of  cloth  airflow,  n. 

O  -  Fmax  ®  1  =  90° 

'  F  @T  -  -90"  1121 

max  ^ 

As  the  altitude  continues  to  increase,  the  force  ratio  R 
approaches  one  indicating  that  at  some  altitude  the  opening  shock 
force  is  essentially  independent  of  the  trajectory  deployment 
angle.  Note  that  the  effect  of  the  value  of  n  is  converging. 
Figure. 6  shows  that  at  high  altitudes  the  inflation  reference  time 
is  basically  no  longer  a  function  of  n  or  k.  The  air  density  is  a 
factor  affecting  the  rate  of  flow  through  the  canopy  surface.  As 
the  density  approaches  zero  the  mass  flow  ratio,  M' ,  approaches 
zero.  See  Section  VII  of  Appendix  A  for  a  discussion  of  Mass  Flow 
Ratio.  Figures  7  through  11  illustrate  the  effects  of  cloth 
airflow  rate  on  the  inflation  reference  time,  inflation  distance 
and  opening  shock  force  as  a  function  of  altitude  for  the 
parameters  of  example  1. 
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COMPARISON  OF  OPENING  SHOCK  FORCES  CALCULATED  BY  THE  METHODS  OF 
REFERENCE  1  AND  THE  GENERAL  TRAJECTORY  TECHNIQUE 

The  opening  shock  force  comparisons  are  limited  to  a 
trajectory  deployment  angle  of  7  =0  .  The  horizontal  trajectory 
deployment  reference  time  is  calculated  via  equation  (14A)  from 
page  A-7  of  Appendix  A. 


2W 

P9VsCdS0 


gpv0 

CdS0 

r  2  W 

/  \1 12 

e 

^Mo  Aso  k  ^  j 

_  \  / 

-  1 

- 

_ 

(14A) 


In  equation  (14A)  T  was  assumed  to  be  zero  and  n  was  assigned  the 
value  of  0.500. 

The  MIL-C-7020,  type  III  canopy  cloth  airflow  constant  is  k=1.46 
and  the  average  pressure  coefficient  is  taken  to  be  c  '.  =  1.7. 
For  a  sea  level  velocity  of  325.25  fps.  p 

t  =  _ 14  x  200 _ _ T  ki  _  . 

0.002378  x  32.2  x  325.25  x  721  58  |_ 


KI 


32.2  x  0.002378  x  4690.83 

2  x  200 


721.58 


399.53  -  962.1  (1.46) 


^0.002378  x  1.7  j1 12 


tQ  =  0.913  sec. 


The  ballistic  mass  ratio  for  t 

o 


M  « 


2W 

pgvstoc0so 


0.913  sec. 


M  = 


_ 2  x  200 _ 

0.002378  x  32.2  x  325.25  x  0.913  x  721.58 


M  =  0.0244 
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St A  LEVEL 


A.  INFLATION  REFERENCE  TIME  (SEC) 


20  KFT 


B.  INFLATION  REFERENCE  TIME  (SEC) 


50  KFT  TOO  KFT 


C.  INFLATION  REFERENCE  TIME  (SEC)  D.  INFLATION  REFERENCE  TIME  (SEC) 


CLOTH 

k 

n 

MIL-C-7020,  TYPE  III 
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0  63246 

MIL  C-8021,  TYPE  1 

0.99469 

0.57403 

MIL-C-8021 ,  TYPE  II 

0.441 17 

0.60492 

FIGURE  6.  EFFECT  OF  CANOPY  CLOTH  CONSTANTS  "k”  and  "n"  ON  THE  PARACHUTE  INFLATION 

REFERENCE  TIME  "T  " 

o 
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INFLATION  REFERENCE  TIME  (SEC) 

FIGURE  7.  EFFECT  OF  ALTITUDE  ON  THE  INFLATION  REFERENCE  TIME  "T  " 
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INFLATION  DISTANCE 

FIGURE  8.  EFFECT  OF  ALTITUDE  ON  THE  INFLATION  DISTANCE  (FEET) 
n  =  0.500, 7  =  0°  t  =  0 
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0  100  200  300  400 


INFLATION  DISTANCE 


FIGURE  9.  EFFECT  OF  ALTITUDE  ON  THE  INFLATION  DISTANCE  (FEET) 

n  =0.632,  y  =  0°  ,  T  =  0 
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FIGURE  10.  VARIATION  OF  STEADY-STATE  DRAG  AND  MAXIMUM  OPENING  SHOCK  FORCE  WITH 
ALTITUDE  FOR  n  =  0.500,  y  =  0°  r  =  0 
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FIGURE  11.  VARIATION  OF  STEADY-STATE  DRAG  AND  MAXIMUM  OPENING  SHOCK  FORCE  WITH 
ALTITUDE  FOR  n  =  0.632,  y  =  Of  r  =  0 
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Since  the  ballistic  mass  ratio  is  less  than  the  limiting 
value  of  =  0.1907,  the  inflation  performance  is  in  the  finite 
mass  range  of  operation  and  t  =  tf.  The  following  formulae  were 
obtained  from  Table  12  of  Reference  1  for  values  of  j=6  and  T =0 . 

The  time  of  occurrence  of  the  maximum  shock  force  during 
canopy  inflation: 


0.7454 


The  maximum  shock  factor: 

y  =  1_6  /  21 M  \6/7 

imax  49  V  4  / 


=  1§ 
'max  aq 


X,  =  0.0560 


'max 


The  steady-state  drag  force: 


fs  =  ypVs2cds0 


Fs  =  -J  x  0.002378  (340)2  721.6 


Fs  =  90761.4  lb. 


The  maximum  shock  force  in  the  horizontal  mode  of  deployment: 


f, 


max 


F max  =  0.0560  x  90761.4 


max 


Fmax  =  5082.6 


max 
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TABLE  5.  COMPARISON  OF  OPENING  SHOCK  FORCES  CALCULATED  FROM 
THE  METHOD  OF  REFERENCE  1  AND  THE  GENERAL  TRAJECTORY 
TECHNIQUE  y  =  0°,  n  =  0.500,  T=  0 


METHOD  OF  REFERENCE  1 

RATIO 

ATL 

Kft 

to 

sec. 

M 

f/to  @  *imax 

Y. 

'max 

Fs 

lb. 

Fmax 

lb. 

GTT 

REF.  1 

> 

K* 

O 

0 

0.914 

0.024 

0.745 

0.0560 

90,761.4 

5,082.6 

1.0078 

o 

—1 

LLt 

> 

H 

2 

3 

0.808 

0.030 

0.768 

0.0672 

83,058.4 

5,581.5 

0.9975 

< 

f— 

CO 

Z 

O 

o 

20 

0.485 

0.086 

0.893 

0.1655 

48,390.3 

8,008.6 

0.9985 

0 

0.914 

0.024 

0.745 

0.0560 

90,761.4 

5,082.6 

1.0048 

JT 

H 

Z 

< 

K 

3 

0.773 

0.030 

0.768 

0.0672 

90,761.4 

6,099.2 

0.9977 

CO 

z 

o 

o 

20 

0.354 

0.086 

0.893 

0.1655 

90,761.4 

15,021.0 

0.9930 

GENERAL  TRAJECTORY  TECHNIQUE  (GTT) 
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A  series  of  opening  shock  force  calculations  were  conducted  for 
constant  velocity  and  constant  dynamic  pressure  profiles  for 
trajectory  deployment  angles  of  7  =0  at  sea  level  and  altitudes 
of  3,000  feet  and  20,000  feet.  The  results  are  summarized  in 
Table  5.  An  examination  of  Table  5  indicates  that  the  forces 
calculated  by  the  method  of  Reference  1  and  the  general  trajectorv 
technique  agree  within  less  than  +  1%,  and  therefore  the  methods 
are  consistent  and  interchangeable. 


DISCUSSION  OF  SYSTEM  PERFORMANCE 

The  listed  performance  of  method  1  shows  the  same  values  for 
M,  t/t  @  x±  max,  and  x.  for  constant  velocity  and  constant 
dynamic  pressure  at  a  designated  altitude.  Yet  the  inflation 
reference  times  and  opening  shock  forces  vary.  How  this  occurs  is 
as  follows: 

Two  major  design  considerations  for  any  parachute  system  are 
the  required  rate  of  descent  (V  )  and  the  altitude  at  which  it  is 
to  be  achieved.  e 


w 


Once  these  decisions  are  firm  the  system,  W/CnSn  has  been 
established.  The  particular  system  weight  to Dbe  recovered  defines 
the  parachute  aerodynamic  size  (CDSn).  The  geometric  size  (Dn) 
and  inflated  shape  and  volume  are  a  result  of  the  type  of  0 
parachute  and  canopy  rate  of  airflow  selected. 

Multiplying  both  sides  of  equation  (14A)  by  V  shows  that  the 
inflation  distance  is  dependent  on  the  parachute  glometry,  cloth 
airflow  properties,  W/CDSo,  and  altitude. 


gpvc 

cns„ 

V  t  =  2W 

-  2W 

e 

A  Mo  -  Aso  k  ( 

,crr  1 

- 1 

s  o 

p g  cdso 

When  the  parachute  type,  canopy  cloth,  and  suspension  line 
length  have  been  selected,  the  only  variable  undefined  is  the 
operational  altitude.  The  inflation  distance  does  not  depend  on 
the  operational  velocity,  but  the  inflation  reference  time  does. 
Rather,  the  inflation  distance  has  been  preordained  by  the  system 
requirements,  size,  and  operational  altitude. 
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The  inflation  distance  is  used  to  determine  the  BMR  which  is 
also  not  a  function  of  velocity,  but  does  depend  on  the 
operational  altitude. 


M 


2W 


PgVstoCDSo 


The  time  of  occurrence  of  the  maximum  shock  factor  during 
the  inflation  process,  and  the  maximum  shock  factor  are  defined  by 
the  BMR. 


1/7 


16_  /21M\  6/7 


The  conclusion  to  be  drawn  from  this  is  that  the  entire 
performance  of  the  parachute  system  was  defined  by  the  initial 
decisions  of  rate  of  descent,  altitude  at  which  the  rate  of 
descent  is  required,  system  weight,  parachute  type,  canopy  cloth, 
and  parachute  geometry  including  suspension  line  length.  The 
performance  varies  with  operational  altitude. 

The  BMR  is  the  scale  parameter  which  determines  what 
percentage  of  the  steady-state  drag  force,  F  ,  is  to  be  realized 
as  opening  shock  force.  As  the  operational  Sensity  altitude  is 
raised  the  inflation  distance  is  reduced  and  the  resulting  BMR  is 
raised.  Selecting  a  constant  velocity-altitude  profile  determines 
the  locus  of  the  variation  of  the  steady-state  drag  force.  A 
constant  dynamic  pressure  locus  of  variation  is  another  convenient 
way  of  determining  how  the  required  line  stretch  velocity  varies 
with  altitude  for  the  purpose  of  providing  a  basis  for  analysis 
and  comparison.  However,  the  percentage  of  the  steady-state  drag 
destined  to  become  opening  shock  force  is  the  same  for  both 
profiles  for  a  given  altitude. 


EXAMPLE  2  -  Compare  the  opening  shock  forces  of  a  ribbon 
parachute  retardation  system  using  the  methods  of  Reference  1  and 
the  general  trajectory  technique. 


a . 

Altitude  =  3,000  ft.  ;  p  = 

0.0021752  slugs  per 

b. 

System  weight  =  4,000  lb. 

c . 

Parachute  diameter  =  35  ft 

;  S  =962.1  ft2 
o 

d. 

Trajectory  velocity  =  340 

fps 

e . 

Trajectory  angle  7 =  0° 

f . 

Initial  drag  area  T =  0 

g- 

j  =  l 
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Parachute  inflation  time: 

nD_ 

r  = 


V. 


14x35  where  n=14  for  ribbon  parachutes 


340 


to  =  t441  sec. 


The  parachute  drag  area: 


CDS0  =  0.5  x962.1 


CDS0  =  481  ft* 


The  ballistic  mass  ratio: 

2W 


M  = 


pgVstoCDSc 


M  = 


2  x  4000 


0.0021753  x  32.2  x  340  x  1.441  x  481.0 


M  =  0.4847 


The  time  of  occurrence  in  the  deployment  process: 


max 


\1 6/  \  3  J 


% 


max  \1 6, 


9\  /2  x  0.4847 


=  0.3197 


The  maximum  shock  factor: 


1  ^  /2M 

—  @Xj  =  ( — 
tQ  max  \ 3 


—  @X; 


max 


v  % 

2  x  0.4847\ 


—  @  X;  =  0.5684 
t  max 

o 
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The  steady  state  drag  force: 

FS  =  J  PVS2  CDS0 

1  9  * 

F  =  —  x  0.0021753  (340  r  481.0 
s  2 

Fs  =  60,474.5  1b. 

The  maximum  shock  force: 

F  =  v.  F 

max  X,max  s 

Fmax  =  0.3197  x60,474.5 


Fmax  =  19.333.7  1b. 


The  computer  solution  of  example  2  is  presented  in  Table 
comparison  of  the  results  with  the  method  of  Reference  1 
in  Table  7. 


6  and  a 
is  shown 
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TABLE  6.  GENERAL  TRAJECTORY  TECHNIQUE  MODE  2  COMPUTER  SOLUTION  OF  EXAMPLE  2 


DEPLOYMENT  CONDITIONS 


SYSTEM  PARAMETERS 


ALTITUDE 
DENSITY 
VELOCITY  VS 
J 

TAU 

ALPHAO 


3000  FT 

0.0021  752  SLUGS/FT3 
340  0  FT/SEC 
1 

0.000 
0.00  DEG 


STEADY-STATE  DRAG  AREA  CDSO  481.00  FT2 

WEIGHT  W  4000.0  LBS 

INFLATION  TIME  TO  1  441  SEC 
STEADY-STATE  DRAG  FORCE  FS  60474.5  LBS 
M  -  0.48460 


TIME  SEC 

VELOCITY 

FT/SEC 

ALPHA  DEG 

RANGE  FT 

S  FT 

DISTANCE 

FT 

SHOCK 

FACTOR 

F  LBS 

0.000 

340.00 

0.00 

0  0 

0.00 

0.00 

0.000E  +  00 

0.0 

0.029 

339.86 

0.16 

9.8 

0.01 

9.80 

1.998E  02 

1208.5 

0.058 

339.44 

0.31 

19.6 

0.05 

19.59 

3.987E  02 

241 1.1 

0.086 

338.75 

0.47 

29.4 

0.12 

29.37 

5.956E  02 

3601,9 

0.115 

337.79 

0.63 

39.1 

0.21 

39.13 

7.896E  02 

4775.2 

0.144 

336.56 

0.79 

48.9 

0.33 

48.85 

9.798E  02 

5925.6 

0.173 

335.06 

0.94 

58.5 

0.48 

58.54 

1.165E  01 

7047.7 

0.202 

333.32 

1.10 

68.2 

0.65 

68.18 

1.345E  01 

8136.8 

0.231 

331.32 

1 .26 

77.8 

0.85 

77.76 

1.519E-01 

9188.4 

0.259 

329.10 

1  42 

87  3 

1.07 

87.28 

1.686E  01 

10198.4 

0.288 

326.64 

1  59 

96.7 

1  32 

96.73 

1.846E  01 

11  163.3 

0.317 

323  97 

1.75 

106.1 

1.59 

106  10 

1.997E  01 

12079.7 

0.346 

321.10 

1.91 

1 15.4 

1.89 

115.39 

2.141E-01 

12945.2 

0.375 

318.04 

2.08 

124.6 

2.21 

124.60 

2.275E  01 

13757.6 

0.404 

314.79 

2.25 

133.7 

2  56 

133.71 

2.400E  01 

14515.2 

0.432 

311  38 

2.42 

142.7 

2.92 

142.73 

2.515E  01 

15216.8 

0.461 

307.82 

2  59 

151.6 

3.31 

151.65 

2.623E  01 

15862.0 

0.490 

304.12 

2  76 

160.4 

3.72 

160.47 

2.720E  01 

16450.5 

0.519 

300.29 

2.94 

169.1 

4.16 

169.18 

2.808E  01 

16982.2 

0.548 

296.34 

3.12 

177.7 

4.61 

177.77 

2.887E  01 

17457.5 

0,576 

292.30 

3.30 

186.2 

5.09 

186.25 

2.956E  01 

17878.0 

0.605 

288.16 

3.48 

194  5 

5.58 

194.61 

3.017E  01 

18244.6 

0.634 

283.95 

3.67 

202.7 

6.09 

202.85 

3.069E  01 

18559.3 

0.663 

279.68 

3.85 

210.8 

6,63 

210.97 

3.113E  01 

18823.2 

0.692 

275.35 

4.05 

218.8 

7.18 

218.97 

3.148E  01 

19038.0 

0.721 

270.98 

4.24 

226.7 

7  74 

226  84 

3.176E  01 

19207.0 

0  749 

266.59 

4.44 

234  4 

8.33 

234.58 

3.197E  01 

19332.6 

0.778 

262.17 

4.64 

242.0 

8.93 

242.20 

3.21 1  E  01 

19417.2 

0.807 

257.72 

4.84 

249.4 

9.55 

249.69 

3.218E  01 

19458.4 

0.823 

255.23 

4.96 

253.6 

9.91 

253  84 

3.219E  01 

19467.5  MAX  FORCE 

0.836 

253.30 

5.05 

256.8 

10.19 

257.05 

3.219E  01 

19467.3 

0.865 

248.88 

5.26 

264.0 

10.84 

264.28 

3.215E  01 

19442.4 

0.894 

244.47 

5.47 

271.0 

1  1.50 

271.39 

3.205E  01 

19384.1 

0.922 

240.08 

5.69 

278.0 

12.18 

278.37 

3.191E  01 

19297.0 

0.951 

235.71 

5.91 

284.8 

12.87 

285.22 

3.172E  01 

19182.9 

0.980 

231.38 

6  14 

291.5 

13  58 

291.95 

3.149E  01 

19044.3 

1 .009 

227.08 

6.37 

298.1 

14.30 

298.55 

3.123E  01 

18883.3 

1 .038 

222.83 

6.61 

304.5 

15  03 

305.03 

3.093E  01 

18702.4 

1.066 

218.63 

6.85 

310.8 

15.77 

31  1.39 

3.060E  01 

18503.7 

1.095 

214.48 

7.09 

317.0 

16.53 

317.63 

3.024E  01 

18289.2 

1.124 

210.38 

7.34 

323.1 

17.30 

323.75 

2.986E  -  0 1 

18060.4 

1 .153 

206.34 

7.59 

329.0 

18.08 

329.76 

2.947E  01 

17819.3 

1 .182 

202.37 

7.85 

334.9 

18.87 

335.64 

2.905E  01 

17567.8 

1 .211 

198.46 

8.11 

340.6 

19  67 

341.42 

2.862E  01 

17307.4 

1.239 

194.61 

8.38 

346.2 

20.48 

347.08 

2.818E  01 

17039.4 

1 .268 

190.83 

8.65 

351.7 

21.31 

352.63 

2.772E  01 

16765.1 

1.297 

187.12 

8.93 

357  1 

22.14 

358  07 

2.726E  01 

16486.1 

1 .326 

183.48 

9.21 

362.3 

22.98 

363.41 

2.679E  01 

16203.1 

1 .355 

179.92 

9.50 

367.5 

23.83 

368.65 

2.632E  01 

15917.6 

1  .384 

176.42 

9.79 

372.6 

24  69 

373.78 

2.585E  01 

15630.5 

1.412 

172.99 

10.09 

377.5 

25.56 

378.81 

2.537E  -  01 

15342.4 

1  441 

169.64 

10.40 

382.4 

26.44 

383.75 

2.489E  01 

15054.0 
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TABLE  7.  COMPARISON  OF  THE  RESULTS  OF  EXAMPLE  2 


tA0  @  Xj 

max 

^max 

METHOD  OF  REFERENCE  1 

0.5684 

19,333.7 

GENERAL  TRAJECTORY 

TECHNIQUE 

0.5711 

19,467.5 

The  variation  of  the  particular  values  determined  by  the  two 
methods  is  less  than  one  percent. 
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CONCLUSIONS 


1 .  The  comparison  of  the  maximum  opening  shock  forces  computed  by 
the  methods  of  reference  1  and  this  report  for  trajectory 
deployment  angles  of  zero  degrees  are  in  agreement  by  less  than 
plus  or  minus  one  percent.  The  two  methods  are  therefore 
consistent  and  interchangeable. 

2.  For  a  given  system  deployed  under  similar  conditions  of 
altitude  and  velocity  the  maximum  opening  shock  force  depends  on 
the  trajectory  angle  at  parachute  inflation.  The  minimum  shock 
force  occurs  in  the  vertical,  away-from-the-earth  trajectory,  and 
the  maximum  shock  force  occurs  in  the  vertical,  toward-the-earth 
trajectory.  As  altitude  increases  the  ratio  of  the  maximum 
opening  shock  force  in  the  vertical  toward-the-earth  trajectory  to 
the  maximum  opening  shock  force  in  the  vertical  away- from- the 
earth  trajectory  approaches  one.  The  opening  shock  force  at  high 
altitudes  is  not  as  dependent  on  trajectory  deployment  angle  as  it 
is  at  low  altitudes. 

3.  As  the  trajectory  deployment  angle  varies  from  a  climb  to  a 
dive  attitude  the  maximum  shock  force  varies  from  minimum  to 
maximum. 

4.  At  low  altitudes  the  cloth  airflow  parameters  k  and  n  affect 
the  canopy  inflation  characteristics.  Extending  the  inflation 
distance  results  in  a  lowering  of  the  maximum  shock  force.  As  the 
altitude  increases,  the  effects  of  the  variation  of  k  and  n  on  the 
inf ldlion  distance  are  minimized  due  to  the  air  density  reduction. 

5.  As  the  altitude  increases  the  maximum  shock  force  increases 
due  to  an  increase  in  the  BMR. 

6.  The  weight-to-drag-area  ratio  is  usually  determined  from  the 
requirement  for  a  given  rate  of  descent  at  a  given  altitude.  This 
ratio,  the  choice  of  parachute  type,  canopy  cloth,  and  suspension 
lir^®  length  define  the  inflation  distance  of  the  system  for  a 
given  altitude.  The  inflation  distance  decreases  as  altitude 
increases.  The  inflation  distance  is  independent  of  deployment 
velocity;  however,  the  inflation  reference  time  is  dependent  upon 
deployment  velocity. 

7.  The  BMR,  which  is  dependent  on  the  inflation  distance  and 
altitude  is  also  independent  of  deployment  velocity.  The  time  of 
occurrence  of  the  maximum  shock  force  during  deployment  and  the 
magnitude  of  the  maximum  shock  factor  are  determined  by  the  BMR 
type  of  parachute  (j),  and  the  initial  drag  area  (  7*  )  .  The 

P ^ oyirie nt  performance  is  independent  of  deployment  velocity. 
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8.  The  BMR  determines  what  percentage  of  the  steady-state  drag 
force  is  to  be  realized  as  opening  shock  force  in  the  particular 
deployment  condition.  Selection  of  a  constant  velocity  or  a 
constant  dynamic  pressure  versus  altitude  profile  determines  the 
magnitude  of  the  steady-state  drag  force  which  determines  the 
opening  shock  force. 

9 .  The  program  mode  1  method  for  solid-cloth  parachutes  can  be 
expanded  to  other  types  of  canopies  when  methods  of  describing  the 
particular  canopy  outflow  have  been  developed. 
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NOMENCLATURE 


Symbols  required  for  this  report  are  listed  on  page  A-16  of 
Appendix  A.  Additional  symbols  not  listed  in  Appendix  A  are  as 
follows : 

j  Time  ratio  exponent  used  in  the  Pflanz  and  modified 

Pflanz  drag  area  equations 

N  Canopy  depth  of  the  inflated  parachute 

R  Ratio  of  the  maximum  opening  shock  force  in  the 

toward-the-earth  trajectory  to  the  opening  shock 
force  in  the  away-f rom-the-earth  trajectory 

ve  Vertical  rate  descent,  7  =  90  degrees,  of  the 

parachute  system  at  the  design  altitude 

7  Trajectory  angle  to  the  horizontal  at  the  time  of 

parachute  suspension  line  stretch.  Loft 
trajectories  are  noted  as  negative  angles.  Dive 
angles  are  noted  as  positive  angles 

T  Ratio  of  parachute  drag  area  at  t=0  to  the  steady- 

state  drag  area 
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A  TECHNIQUE  FOR  THE  CALCULATION  C?  THE  CPENIUG-SHOCK  FORCES  FOR 
SEVERAL  TYPES  OF  SOLIL  CLOTH  PARACHUTES 

W.  P.  Ludtfce 

Naval  Ordnance  Laboratory 
Silver  Spring,  Maryland 


Abstract 


An  analytical  method  of  calculating 
parachute  opening-shock  forces  based  upon 
wind-tunnel  derived  drag  area  time  signa¬ 
tures  of  several  solid  cloth  parachute 
types  in  conjunction  with  a  scale  factor 
and  retardation  system  steady-state  param¬ 
eters  has  been  developed.  Methods  of 
analyzing  the  inflation  time,  geometry, 
cloth  airflow  properties  and  materials 
elasticity  are  included.  The  effects  of 
mass^  ratio  and  altitude  on  the  magnitude 
and  time  of  occurrence  of  the  maximum 
opening  shock  are  consistent  with  observed 
field  test  phenomena. 

I.  Introduction 


altitude  increases,  effects  of  altitude  on 
opening-shock  force,  finite  and  infinite 
mass  operation,  and  inflation  distance. 

II.  Development  of  Velocity  Ratio  and 
Force  Ratio  Equations  During 
the  Unfolding  Phase  of 
Parachute  Deployment 

The  parachute  deployment  would  take 
place  ir.  a  horizontal  attitude  in  accord¬ 
ance  with  Newton's  second  law  of  motion. 

T,  F  =  ma 


-LAcns  =  * 


dV 
d  t 


In  1965,  the  Naval  Ordnance  Laboratory 
(NOL)  was  engaged  in  a  project  which 
utilized  a  35-foot-diameter,  10-percent 
extended-skirt  parachute  (type  T-10)  as 
the  second  stage  of  a  retardation  system 
for  a  250-pound  payload.  Deployment  of 
the  T-10  parachute  was  to  be  accomplished 
at  an  altitude  of  100,000  feet.  In  this 
rarefied  atmosphere,  the  problem  was  to 
determine  the  second  stage  deployment 
conditions  for  successful  operation.  A 
search  of  available  field  test  information 
indicated  a  lack  of  data  on  the  use  of 
solid  cloth  parachutes  at  altitudes  above 
30,000  feet. 

The  approach  to  this  problem  was  as 
follows:  Utilizing  existing  wind-tunnel 

data,  low-altitude  field  test  data,  and 
reasonable  assumptions,  a  unique  engi¬ 
neering  approach  to  the  inflation  time 
and  opening-shock  problem  was  evolved  that 
provided  satisfactory  results.  Basically, 
the  method  combines  a  wind-tunnel  derived 
drag  area  ratio  signature  as  a  function 
of  deployment  time  w ith  a  scale  factor  and 
Newton's  second  law  of  motion  to  analyze 
the  velocity  and  force  profiles  during 
deployment.  The  parachute  deployment 
sequence  is  divided  into  two  phases.  The 
first  phase,  called  "unfolding  phase," 
where  the  canopy  is  undergoing  changes  in 
shape,  is  considered  to  be  inelastic  as 
the  parachute  inflates  initially  to  its 
steady-state  aerodynamic  size  for  the 
first  time.  At  this  point,  the  "elastic 
phase"  is  entered  where  it  is  considered 
that  the  elasticity  of  the  parachute 
materials  enters  the  problem  and  resists 
the  applied  forces  until  the  canopy  has 
reached  full  inflation. 

The  developed  equations  are  in  agree¬ 
ment  with  the  observed  performance  of 
solid  cloth  parachutes  in  the  field,  such 
as  the  decrease  of  inflation  time  as 


It  was  recognized  that  other  factors,  such 
as  included  air  mass,  apparent  mass,  and 
their  derivatives,  also  contribute  forces 
act]'  z  or-  the  system.  Since  definition 
of  txxese  parameters  was  difficult,  the 
analysis  was  conducted  in  the  simplified 
form,  shown  above.  Comparison  of  calculated 
results  and  test  results  indicated  that 
the  omitted  terms  have  a  small  effect. 


t  V 


Multiplying  the  right-hand  side  of  equa¬ 
tion  (1)  by 


1  = 


Yt  C.SA 
s  o  D  o 

s  o  D  o 


and  rearranging 


i_  r  cps 

jo  Vo 


-2W 


PgV  t  C-.S 
S  O  D  o 


(2) 


In  order  to  integrate  the  left-hand  term 
of  equation  (2),  the  drag  area  ratio  must 
be  defined  for  the  type  of7 parachute  under 
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analysis  as  a  function  of  deployment 

v.i-’aTT.'-.p/i  f  i  r-  r.  *f- 

-  "  ••-  '  -l~<-  - j  • 


TIME  RATIO  t/tQ 

FIG.  1  TYPICAL  INFINITE  MASS  FORCE-TIME  HISTORY  OF  A 
SOLID  CLOTH  PARACHUTE  IN  A  WIND  TUNNEL 


{LB) 

120 


0.10  0.20  0.30 

TIME  (SEC) 

REPRODUCED  FROM  REFERENCE  (1) 

FIG.  3  TYPICAL  FORCE-TIME  CURVE  FOR  A  10%  EXTENDED 

SKIRT  PARACHUTE  UNDER  INFINITE  MASS  CONDITIONS. 
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TIME  (SEC) 

REPRODUCED  FROM  REFERENCE  (1) 
FIG.  4  TYPICAL  FORCE-TIME  CURVE  FOR  A  PERSONNEL 
GUIDE  SURFACE  PARACHUTE  UNDER  INFINITE 
MASS  CONDITIONS 


Figure  1  illustrates  a  typical  solid 
cloth  parachute  wind-tunnel  infinite  mass 
force-time  history  after  snatch.  In 
infinite  mass  deployment*  the  maximum  size 
ar.d  maximum  shock  force  occur  at  the  time 
of  full  inflation*  tf.  However *  tf  is 
inappropriate  for  analysis  since  it  is 
dependent  upon  the  applied  load*  structural 
strength*  and  materials  elasticity.  The 
reference  time*  tc*  where  the  parachute  has 
attained  its  steady-state  aerodynamic  size 
for  the  first  time,  is  used  as  the  basis  for 
performance  calculations. 

At  any  instant  during  the  unfolding 
phase*  the  force  ratio  F/Fs  can  be  deter¬ 
mined  as  a  function  of  the  time  ratio, 
t/.c. 


~  s 


D~o 


Since  the  wind-tunnel  velocity  and  density 
are  constant  during  infinite  mass  deployment 


CDS 

CDSo 


FIG.  2  TYPICAL  FORCE-TIME  CURVE  FOR  A  SOLID  FLAT  PARA¬ 
CHUTE  UNDER  INFINITE  MASS  CONDITIONS. 


FIG.  5  TYPICAL  FORCE-TIME  SIGNATURE  FOR  THE  ELLIPTICAL 
PARACHUTE  UNDER  INFINITE  MASS  CONDITIONS 


TI^C  fSECI 


FIG.  6  TYPICAL  FORCE-TIME  SIGNATURE  FOR  THE  RING 

SLOT  PARACHUTE  20%  GEOMETRIC  POROSITY  UNDER 
INFINITE  MASS  CONDITIONS 


Infinite  mass  opening- shock  signatures 
cf  several  types  of  parachutes  are  pre¬ 
sented  in  figures  2  through  6.  Analysis 
of  chese  signatures  using  the  force  ratio* 
F/?s*  -  time  ratio*  t/tQ*  technique 
indicated  a  similarity  in  the  performance 
of  the  various  solid  cloth  types  of 
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parachutes  which  were  examined.  The  geo¬ 
metrically  porous  ring  slot  parachute 
displayed  a  completely  different  signature, 
as  was  expected.  These  data  are  illus¬ 
trated  in  Figure  7-  If  an  initial  boundary 


CD_S 

cDs0 


+  277(1 


+  r/2 


(5) 


FIG.  7  DRAG  AREA  RATIO  VS.  TIME  RATIO 

condition  of  C^S/C^S  -  0  at  time  t/t  -  0 
D  Do  '  o 

is  assumed,  then,  the  data  can  be  approxi¬ 
mated  by  fitting  a  curve  of  the  form 


At  the  time  that  equation  (5)  was  ascer¬ 
tained,  it  suggested  that  the  geometry  of 
the  deploying  parachute  was  independent 
of  density  and  velocity.  It  v:as  also 
postulated  that  although  this  expression 
had  been  determined  for  the  infinite  mass 
condition,  it  would  also  be  true  for 
the  finite  mass  case.  This  phenomenon 
has  since  been  independently  observed  and 
confirmed  oy  Berndt  and  Be  Teese  in 
reference  (£). 

Since  the  drag  area  ratio  was  determined 
from  actual  parachute  deployments,  it  v/as 
assumed  that  the  effects  of  apparent  mass 
and  included  mass  on  the  deployment  force 
history  were  accommodated. 

The  right-hand  term  of  equation  (2) 
contains  the  expression 


PSVst0CDS0 

This  term  can  be  visualized  as  shown  in 
Figure  8  to  be  a  ratio  of  the  retarded 
mass  (including  the  parachute)  to  an 
associated  mass  of  atmosphere  contained 
in  a  right  circular  cylinder  which  is 
generated  by  moving  an  inflated  parachute 
of  area  C^SQ  for  a  distance  equal  to  the 
product  of  Vst0  through  an  atmosphere  of 
density,  o. 


A  more  realistic  drag  area  ratio  expres¬ 
sion  was  determined  which  includes  the 
effect  of  initial  area  at  line  stretch. 


CDS 


CDSo 


where  r\  is  the  ratio  of  the  projected 
mouth  area  at  line  stretch  to  the  steady- 
state  projected  frontal  area.  Expanding 
equation  (4) 


FIG.  8  VISUALIZATION  OF  THE  MASS  RATIO  CONCEPT 


The  mass  ratio,  M,  is  the  scale  factor 
which  controls  the  velocity  and  force 
profiles  during  parachute  deployment. 
Substituting  K  and  Cp3/CpSo  into  equa¬ 
tion  (2),  integrating,  and  solving  for 

Y/Y 
/  s 


V_ 

V5 


1 


1  ( 
M 


+  *?2 


t 

to 


(7) 
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The  instantaneous  shock  factor  is  defined 
as 


x.  = 

i 


F 


Figures  9  and  10  illustrate  the  velocity 
and  force  profiles  generated  from  equations 
(7)  and  (8)  for  initial  projected  area 
ratios  of  t)  =  0,  and  0.2  with  various  mass 
ratios . 


IV.  Methods  for  Calculation  of 
the  Reference  Time.  tQ 


If  the  altitude  variation  during  deployment 
is  small,  then,  the  density  may  be  consid¬ 
ered  as  constant 


xi  = 


CDS 

CDSo 


(r) 


from  equations  (5)  and  (7) 


(1  -  T})2 


1+! 

Ri -n)2/ tV 

M 

?  u 

+  2r?(1  -V)[~J  +1?* 


(8) 


III.  Maximum  Shock  Force  and  Time  of 
Occurrence  During  the  Unfolding  Phase 

The  time  of  occurrence  of  the  maximum 
instantaneous  shock  factor,  Xj_,  is  diffi¬ 
cult  to  determine  for  the  general  case. 
However,  for  ^  =  0,  the  maximum  shock 
factor  and  time  of  occurrence  are  readily 
calculated.  For  r\  =  0 


The  ratio  concept  is  an  ideal  method  to 
analyze  the  effects  of  the  various  param¬ 
eters  on  the  velocity  and  force  profiles 
of  the  opening  parachutes;  however,  a 
means  of  calculating  tQ  is  required  before 
specific  values  can  be  computed.  Methods 
for  computing  the  varying  mass  flov;  into 
the  inflating  canopy  mouth,  the  varying 
mass  flov;  out  through  the  varying  inflated 
canopy  surface  area,  and  the  volume  ex’ 
air,  which  must  be  collected  during 
the  inflation  process  are  required. 

Figure  11  represents  a  solid  cloth- 
type  parachute  canopy  at  some  instant 
during  inflation.  At  any  given  instant, 
the  parachute  drag  area  is  proportional 
to  the  maximum  inflated  diameter.  Also, 
the  maximum  diameter  in  conjunction  with 
the  suspension  lines  determines  the  inflow 
mouth  area  (A- A)  and  the  pressurized 
canopy  area  (B-B-B).  This  observation 
provided  the  basis  for  the  following 
assumptions.  The  actual  canopy  shape  Is 
cf  minor  importance. 

a.  The  ratio  of  the  instantaneous 
mouth  inlet  area  to  the  steady- state 
mouth  area  is  in  the  same  ratio  as  the 
instantaneous  drag  area. 


Setting  the  derivative  of  xi  with  respect 

to  time  equal  to  zero  and  solving  for 

t/t  at  x. 

o  i  max 


1 


and  the  maximum  shock  factor  is 

*  « -  f;  iff  <”> 


_\l__  CDS 

^Mo  CBSo 

b.  The  ratio  of  the  instantaneous 
pressurized  cloth  surface  area  to  the 
canopy  surface  area  is  in  the  same  ratio 
as  the  instantaneous  drag  area. 


c.  Since  the  suspension  lines  in  the 
unpressurized  area  of  the  canopy  are 
straight,  a  pressure  differential  has  not 
developed,  and,  therefore,  the  net  air¬ 
flow  in  this  zone  is  zero. 

Based  on  the  foregoing  assumptions, 
the  mass  flov;  equation  can  be  written 


Equations  (9)  and 
of  M  *  A  (0.19), 

of  H,  the  maximum 
the  elastic  phase 


(10)  are  valid  for  valuta 
since  for  larger  values 

shock  force  occurs  in 
of  inflation. 


dm  =  m  in  flov;  -  m  out  flov; 


p  « =  pVA> 


-  PA sp 
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FIG.  9  EFFECT  OF  INITIAL  AREA  AND  MASS  RATIO  ON  THE 
SHOCK  FACTOR  AND  VELOCITY  RATIO  DURING  THE 


UNFOLDING  PHASE  FOR 
M-CD 


M  =4 
M  =  2 

M  =  1 
M=  0.8 

M  =  0.6 

M  =  0.4 


M=  0.2 


M  =  0.1 
M  «  0.08 
M  =  0.06 
M  =  0.04 

M=  0.02 


t/tQ  UNFOLDING  TIME  RATIO  t/t 

FIG.  10  EFFECT  OF  INITIAL  AREA  AND  MASS  RATIO  ON  THE 
SHOCK  FACTOR  AND  VELOCITY  RATIO  DURING  THE 
UNFOLDING  PHASE  FOR  17=0.2 
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From  equation  (3) 


CDS 

CDSo 


■(t)  * 


=  ( —  I  5  for  T|  =  0 


From  equation  (7) 


V  =  - 2 - 

1  +  —  (-A 

7M  \tQ  ) 


;  q  =  o 


From  equation  (2o) 


V 


2n 


Integrating: 


Measured  values  of  n  indicate  a  data 
range  from  0.574  through  6.771.  A 
convenient  solution  to  the  reference  time 
equation  evolves  when  n  is  assigned  a 
value  of  1/2.  Integrating  equation  (13) 
and  using 


V  t  M  = 
s  o 


2W 


gPCDS0 


LET  K, 


gpvc 

2W 


CDSo 


/CpP \ 1/2 

Vo  -  ASok(v 


) 


l4w 


gPVsCDSo 


Kl  n 
e  -  1 


(14) 


Equation  (l4)  expresses  the  unfolding 
reference  time,  tQ,  in  terms  of  mass, 
altitude,  snatch  velocity,  airflow  char¬ 
acteristics  of  the  cloth,  and  the  steady- 
state  parachute  geometry.  Note  that  the 
term  gpVo/W  is  the  ratio  of  the  included 
air  mass  to  the  mass  of  the  retarded 
hardware.  Multiplying  both  sides  of 
equation  (14)  by  Vs  demonstrates  that 


V  t  =  a  constant  which  is 

a  function  of  altitude 


Figures  12  and  13  indicate  the  para¬ 
chute  unfolding  time  and  unfolding 
distance  for  values  of  n  =  1/2  and  n  = 
0.63246.  Note  the  variation  and  conver¬ 
gence  with  rising  altitude.  The  opening- 
shock  force  is  strongly  influenced  by 
the  inflation  time.  Because  of  this,  the 
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value  of  tQ  calculated  by  using  a  real¬ 
istic  value  of  n  should  be  used  in  the 
lower  atmosphere. 

As  an  example  of  this  method  of  opening- 
shock  analysis,  let  us  examine  the  effect 
of  altitude  on  the  opening-shock  force  of 
a  T-10-type  parachute  retarding  a  ^00- 
pound  weight  from  a  snatch  velocity  of 
Vs  =  400  feet  per  second  at  sea  level. 
Conditions  of  constant  velocity  and  con¬ 
stant  dynamic  pressure  are  investigated. 

The  results  are  presented  in  Figure  14. 

At  low  altitudes,  the  opening- shock  force 
is  less  than  the  steady-state  drag  force: 
however,  as  altitude  rises,  the  opening 
shock  eventually  exceeds  the  steady-state 
drag  force  at  some  altitude.  This  trend 
is  in  agreement  with  field  test  observations. 

V.  Correction  of  tQ  for 
Initial  Area  Effects 


The  unfolding  reference  time,  tQ, 
calculated  by  the  previous  methods  assumes 
that  the  parachute  inflates  from  zero  drag 
area.  In  reality,  a  parachute  has  a  drag 


UNFOLDING  TIME  (SEC.) 


FIG.  12  EFFECT  OF  ALTITUDE  ON  THE  UNFOLDING  TIME 
"to"  AT  CONSTANT  VELOCITY  AND  CONSTANT 
DYNAMIC  PRESSURE  FOR  n  =  1/2  AND  n  =  0.63296 


CLOTH  MIL  C-7020,  TYPE  III 

35  FOOT;  10%  EXTENDED  SKIRT  CANOPY 


PIG.  13  EFFECT  OF  ALTITUDE  ON  THE  UNFOLDING  DIS¬ 
TANCE  AT  CONSTANT  VELOCITY  AND  CONSTANT 
DYNAMIC  PRESSURE  FOR  n  -  1/2  AND  n  =  0.63246. 


FIG.  14  VARIATION  OF  STEADY-STATE  DRAG,  F  ,  AND  MAXI¬ 
MUM  OPENING  SHOCK  WITH  ALTITUDE  FOR  CONSTANT 
VELOCITY  AND  CONSTANT  DYNAMIC  PRESSURE 


area  at  the  beginning  of  inflation.  Once 
zQ  has  been  calculated,  a  correction  can 
be  applied,  based  upon  what  is  known  about 
the  initial  conditions. 
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Case  A  -  When  the  initial  projected 
area  is  known 


'■£) 


1/3 


calculated 


than  the  maximum  extensibility,  e 
the  materials.  n 


of 


Development  of  the  analysis  in  the 
elastic  phase  of  inflation  is  similar  to 
the  technique  used  in  the  unfolding  phase. 
Newton* s  second  lav;  of  motion  is  used* 
together  with  the  drag  area  ratio  signa¬ 
ture  and  mass  ratio 


=  (_LY 

CDSo  Vj 


corrected 


calculated 


(15) 


Case  B  -  When  the  initial  drag  area 
is  known 


C  S  / 1  \  ^ 
D  i 


CDSo 


-© 


h 


,C  S  \1/6 

\cDso/ 


calculated 


which  is  still  valid,  as  shown  in  Figure  7 
t  _  V 


1 1  (t)  at=v4 

o 


^d V 


Integrating  and  solving  for  — 

Vo 


Vs  V3 

VD  7M 


(tr 


(17) 


where  —  is  the  velocity  ratio  of  the 


corrected 


/C  s  \1/6' 

- 

VCDSo/ 


unfolding  process  at  time  t  =»  tQ. 


calculated 

(16) 


The  mass  ratio  should  now  be  adjusted  for 
the  corrected  tQ  before  velocity  and  force 
profiles  are  determined. 

VI.  Opening-Shock  Force,  Velocity  Ratio, 
and  Inflation  Time  During  the 
Elastic  Phase  of  Parachute  Inflation 

The  mass  ratio,  M,  is  an  important 
parameter  in  parachute  analysis.  For 
values  of  M  <<4/21,  the  maximum  opening- 
shock  force  occurs  early  in  the  inflation 
process,  and  the  elastic  properties  of 
the  canopy  are  not  significant.  As  the 
mass  ratio  approaches  M  =  4/21,  the 
magnitude  of  the  opening-shock  force 
increases,  and  the  time  of  occurrence 
happens  later  in  the  deployment  sequence. 
For  mass  ratios  M  >  4/21,  the  maximum 
shock  force  will  occur  after  the  reference 
time,  tD.  Parachutes  designed  for  high 
mass  ratio  operation  must  provide  a 
structure  of  sufficient  constructed 
strength,  Fc,  so  that  the  actual  elon¬ 
gation  of  the  canopy  under  load  is  less 


1  +  - 

l1  -  nlf  +  a  U  -  a)  +  2 

M 

7  2 

(18) 


Ihe  instantaneous  shock  factor  in  the 
2lastic  phase  becomes 


x.  =  -°i 

CDSo\Vs/ 

fe)6 


(19) 


The  end  point  of  the  inflation  process 
depends  upon  the  applied 'loads,  elasticity 
of  the  canopy,  and  the  constructed  strength 
of  the  parachute.  A  linear  load  elongation 


i  ♦  -1- 1 

|c+  1 

^  ' 

1 

H 

_ 1  ! 

Vo  71'1  | 

LVto ' 
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relationship  is  utilized  to  determine  the 
maximum  drag  area. 


e 

F 


max 


The  next  force  in  the  series  at  constant  q 

A1 

F1  =  X  F  — 

1  o  s  . 


F  e 


max 


(20) 


wnere 


-  =  (i  +  o2 


The  force,  F,  is  initially  the  instanta¬ 
neous  force  at  the  end  of  the  unfolding 
process 


F  =  X  F 
o  s 


(21) 


Subsequent  elongations  in  the  system  can 
be  shown  to  be 


=  g  (1  +  £  V 
1  o  v  o' 


where  XQ  is  the  shock  factor  of  the  e2  =  eo  +  Go  +  eG)2)2 

unfolding  phase  at  t  =  tD 


Xo  =- 


1  + 


A[ii. 

M  L 


—al-  +  -  a)  +  r,^ 


(22) 


Since  the  inflated  shape  is  defined,  the 
drag  coefficient  is  considered  to  be 
constant,  and  the  instantaneous  force  is 
proportional  to  the  dynamic  pressure  and 
projected  area.  The  maximum  projected 
area  would  be  developed  if  the  dynamic 
pressure  remained  constant  during  the 
elastic  phase.  Under  very  high  mass 
ratios,  this  is  nearly  the  case  over  this 
very  brief  time  period;  but  as  the  mass 
ratio  decreases,  the  velocity  decay  has 
a  more  significant  effect.  The  simplest 
approach  for  all  mass  ratios  is  to 
determine  the  maximum  drag  area  of  the 
canopy  as  if  elastic  inflation  had 
occurred  at  constant  dynamic  pressure. 

Then  utilizing  the  time  ratio  determined 
as  an  end  point,  intermediate  shock 
factors  can  be  calculated  from  equation 
(19)  and  maximum  force  assessed. 

The  initial  force,  X0FS,  causes  the 
canopy  to  increase  in  projected  area.  The 
new  projected  area  in  turn  increases  the 
total  force  on  the  canopy  which  produces 
a  secondary  projected  area  increase.  The 
resulting  series  of  events  are  resisted 
by  the  parachute  materials.  The  parachute 
must,  therefore,  be  constructed  of  suffi¬ 
cient  strength  to  prevent  the  elongation 
of  the  materials  from  exceeding  the 
maximum  elongation. 


G  = 
O 


X  F 
o  s 


(23) 


The  required  canopy  constructed  strength 
can  be  determined  for  a  given  set  of 
deployment  conditions.  The  limiting  value 
of  the  series  (e^)  determines  the  end 
point  time  ratio. 


(!  +  h)2 

(i  ♦ 


{2b) 


Figure  15  illustrates  the  maximum  drag 
area  ratio  as  a  function  of  gq. 


o 

io 


FIG.  15  MAXIMUM  DRAG  AREA  RATIO  VS.  INITIAL 
ELONGATION 
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VII.  Application  of  Cloth  Permeability 
to  the  Calculation  of  the  Inflation 
Time  of  Solid  Cloth  Parachutes 


The  mass  out flow  through  the  pressur¬ 
ized  region  of  an  inflating  solid  cloth 
parachute  at  any  instant  is  dependent  upon 
the  canopy  area  which  is  subjected  to  air¬ 
flow  and  the  rate  of  airflow  through  that 
area.  The  variation  of  pressurized  area  as 
a  function  of  reference  time,  tQ,  was 
earlier  assumed  to  be  proportional  to  the 
instantaneous  drag  area  ratio,  leaving  the 
rate- of- airflow  problem  to  solve.  The 
permeability  parameter  of  cloth  was  a 
natural  choice  for  determining  the  rate  of 
airflow  through  the  cloth  as  a  function  of 
pressure  differential  across  the  cloth. 
Heretofore,  these  data  have  been  more  of 
a  qualitative,  rather  than  quantitative, 
value.  A  new  method  of  analysis  was 
developed  wherein  a  generalized  curve  of 
the  form  P  =  k(AP)n  was  fitted  to  cloth 
permeability  data  for  a  number  of  different 
cloths  and  gives  surprisingly  good  agree¬ 
ment  over  the  pressure  differential  range 
of  available  data.  The  pressure  differ¬ 
ential  was  then  related  to  the  trajectory 
conditions  to  give  a  generalized  expression 
which  can  be  used  in  the  finite  mass  ratio 
range,  as  well  as  the  infinite  mass  case. 
The  permeability  properties  were  trans¬ 
formed  into  a  mass  flow  ratio,  MT ,  which 
shows  agreement  with  the  effective  porosity 
concept . 

Measured  and  calculated  permeability 
pressure  data  for  several  standard  cloths 
are  illustrated  in  Figure  16.  This  method 
has  been  applied  to  various  types  of  cloth 
between  the  extremes  of  a  highly  permeable 
3-momme  silk  to  a  relatively  impervious 
parachute  pack  container  cloth  with 
reasonably  good  results,  see  Figure  17. 

The  canopy  pressure  coefficient,  Ch, 
is  defined  as  the  ratio  of  the  pressure 
differential  across  the  cloth  to  the 
dynamic  pressure  of  the  free  stream. 


AP  =  P(internal)  P(external) 
Q  1/2  pV2 


(25) 

where  V  is  based  on  equation  (7). 

The  permeability  expression,  P  =  k(AP)n 
becomes 


(25) 


Although  some  progress  has  been  made 
by  Melzig  and  others  on  the  measurement  of 
the  variation  of  the  pressure  coefficient 
on  an  actual  inflating  canopy,  this 
dimension  and  its  variation  with  time  are 
still  dark  areas  at  the  time  of  this 
writing.  At  the  present  time,  a  constant 
average  value  of  pressure  coefficient  is 


used  in  these  calculations.  Figure  18 
presents  the  effect  of  pressure  coefficient 
and  altitude  on  the  unfolding  time  for 
constant  deployment  conditions. 

It  is  well  known  that  the  inflation 
time  of  solid  cloth  parachutes  decreases 
as  the  operational  altitude  increases. 

This  effect  can  be  explained  by  consid¬ 
ering  the  ratio  of  the  mass  outflow 
through  a  unit  cloth  area  to  the  mass 
inflow  through  a  unit  mouth  area. 


M1  =  mass  flow  ratio  =  mass 

mass  inflow 


where 

mass  outflow  =  p  (per.  ft2  cloth  area) 

P  ft  -sec 


and 


mass  inflow  =  V 


p  ft2 


slu^s  __  (per  inflow  area) 


REPRODUCED  FROM  REFERENCE  {4) 


FIG.  16  NOMINAL  POROSITY  OF  PARACHUTE  MATERIAL 
VS  DIFFERENTIAL  PRESSURE. 
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Therefore,  the  mass  flow  ratio  becomes 


Effective  porosity,  C,  is  defined  as 
the  ratio  of  the  velocity  through  the 
cloth,  u,  to  a  fictitious  theoretical 
velocity,  v,  which  villi  produce  the  par¬ 
ticular  AP=  l/2ov  . 


o  20  40  60  80  100 


AP  LB/FT2 

FIG.  17  COMPARISON  OF  MEASURED  AND  CALCULATED 

PERMEABILITY  FOR  RELATIVELY  PERMEABLE  AND 
IMPERMEABLE  CLOTHS 


AVERAGE  CANOPY  PRESSUR E  COEFFICI ENT 
DURING  INFLATION  INCLUDING  THE  VENT 


CP 


FIG.  18  EFFECT  OF  PRESSURE  COEFFICIENT  AND  ALTITUDE 
ON  THE  UNFOLDING  TIME. 


effective  porosity,  C  =  —  (28) 


Comparison  of  the  mass  flow  ratio  and 
previously  published  effective  porosity 
data  is  shown  in  Figure  19.  The  effects 
of  altitude  and  velocity  on  the  mass  flow 
ratio  are  presented  in  Figures  20,  and  21 
for  constant  velocity  and  constant 
altitude.  The  decrease  of  cloth  perme¬ 
ability  with  altitude  is  evident. 

The  permeability  constants  nkri  and  !'nM 
can  be  determined  from  the  permeability 
pressure  differential  data  as  obtained 
from  an  instrument  such  as  a  Frazier 
Permeameter .  Two  data  points,  ,:A”  and 


REPRODUCED  FROM  REFERENCE  (4) 


FIG.  19  THE  EFFECTIVE  POROSITY  OF  PARACHUTE 
MATERIALS  VS.  Dl FFERENTIAL  PRESSURE 
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"B/1  are  selected  in  such  a  manner  that 
point  "A"  is  in  a  low-pressure  zone  below 
the  knee  of  the  curve,  and  point  "B"  is 
located  in  the  upper  end  of  the  high- 
pressure  zone,  as  shown  in  Figure  22 . 

The  two  standard  measurements  of  1/2 
inch  of  water  and  20  inches  of  water 
appear  to  be  good  data  points  if  both  are 


FIG.  22  LOCATION  OF  DATA  POINTS  FOR  DETERMINATION 
OF  "k"  AND  "n" 


available  on  the  same  sample.  Substi¬ 
tuting  the  data  from  points  MAIf  and  ,rB” 

into  P  =  k(AP)n  : 


(29) 


k  = 


(APB)n 


(30) 


FIG.  20  EFFECT  OF  ALTITUDE  ON  MASS  FLOW  RATIO  AT 
CONSTANT  VELOCITY 


FIG.  21  EFFECT  OF  VELOCITY  ON  MASS  FLOW  RATIO  AT 
CONSTANT  DENSITY 


VIII.  Determination  of  the 
Parachute  Included  Volume 
and  Associated  Air  Mass 

Before  the  reference  time,  tG,  and 
inflation  time,  tp,  can  be  calculated, 
the  volume  of  atmosphere,  VG  ,  which  is 
to  be  collected  during  the  inflation 
process  must  be  accurately  known.  This 
requirement  dictates  that  a  realistic 
inflated  canopy  shape  and  associated 
volume  of  atmosphere  be  determined. 

Figure  23  was  reproduced  from  reference 
(5).  The  technique  of  using  lampblack 
coated  plates  to  determine  the  airflow 
patterns  around  metal  models  of  inflated 
canopy  shapes  was  used  by  the  investigator 
of  reference  (5)  to  study  the  stability 
characteristics  of  contemporary  parachutes, 
i.e.,  1943.  A  by-product  of  this  study 
is  that  it  Is  clearly  shown  that  the 
volume  of  air  within  the  canopy  bulges 
out  of  the  canopy  mouth  (indicated  by 
arrows )  and  extends  ahead  of  the  canopy 
hem. ^  This  volume  must  be  collected  during 
the  inflation  process.  Another  neglected, 
but  significant,  source  of  canopy  volume 
exists  in  the  billowed  portion  of  the  gore 
panels . 
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HEM  I  SPHERE 


VENT  PARACHUTE 

REPRODUCED  FROM  REFERENCE  (5) 

FIG.  23  AIRFLOW  PATTERNS  SHOWING  AIR  VOLUME 
AHEAD  OF  CANOPY  HEM 


The  steady-state  canopy  shape  has  been 
observed  in  wind-tunnel  and  field  tests  to 
be  elliptical  in  profile.  Studies  of  the 
inflated  shape  and  included  volume  of 
several  parachute  types  (flat  circular., 

10  percent  extended  shirty  elliptical, 
hemispherical,  ring  slot,  ribbon,  and 
cross)  are  documented  in  references  (6) 
and  (7).  These  studies  demonstrated  that 
the  steady-state  profile  shape  of  inflated 
canopies  of  the  various  types  can  be 
approximated  to  be  two  ellipses  of  common 
major  diameter,  2a,  and  dissimilar  minor 
diameters,  b  and  bT,  as  shown  in  Figure  24. 
It  was  also  shown  that  the  volume  of  the 
ellipsoid  of  revolution  formed  by  revolving 
the  profile  shape  about  the  canopy  axis 
was  a  good  approximation  of  the  volume  of 
atmosphere  to  be  collected  during  canopy 
inflation  and  included  the  air  volume 
extended  ahead  of  the  parachute  skirt  hem 
together  with  the  billowed  gore  volume. 


BILLOWED  GORE 
CONTOUR 
A 


EQUIVALENT  BILLOW  VOLUME 
EXTENDED  AIR  VOLUME 


REPRODUCED  FROM  REFERENCE  (6) 


FIG.  24  PARACHUTE  CROSS  SECTION  NOMENCLATURE 


Tables  I  and  II  are  summaries  of  test 
results  reproduced  from  references  (6)  and 
(7) ?  respectively,  for  the  convenience  of 
the  reader. 
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TABLE  I  SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  12-GORE  AMD  16-GORE  CONFIGURATIONS 


Parachute 

No.  of 

Suspension 

Velocity 

Scale  Factor,  K 

N 

Axes  Ratio 

Volume  in 

3 

V 

Type 

Gores 

Line  Length 

2a 

2a 

2a 

2a 

- 

b 

b 

b  b' 

o 

inches 

mph 

fps 

D 

a 

Df 

Dr 

L 

a 

a 

a  5 

Vh 

vc 

V 

o 

VH 

Flat  Circular 

12 

34 

50 

73 

.645 

.650 

.856 

.6115 

.8817 

1.4932 

4476 

4481 

6980 

1.56 

16 

34 

50 

73 

.663 

.669 

.820 

.5558 

.9039 

1.4597 

4450 

4100 

7325 

1.65 

10%  Extended 

12 

34 

100 

147 

.663 

.652 

.881 

.6424 

.8860 

1.5284 

3928 

4400 

6783 

1.73 

Skirt 

16 

34 

17 

25 

.654 

.640 

.785 

.5580 

.8502 

1.4082 

4051 

3920 

6197 

1.53 

Elliptical 

12 

34 

75 

110 

.916 

812 

.5626 

.9657 

1.5283 

3322 

5405 

16 

34 

17 

25 

.875 

.800 

.6169 

8163 

1.4332 

2726 

4405 

Hemispherical 

12 

34 

125 

183 

.996 

1.254 

1.0005 

.9080 

1.9085 

6224 

8666 

16 

34 

75 

110 

994 

1  185 

.9129 

.9380 

1.8509 

5921 

8370 

Ringslot 

12 

34 

25 

37 

.607 

.654 

.853 

.6566 

.8735 

1.530 

3800 

3650 

5903 

1.55 

16%  Geometric 

12 

34 

100 

147 

.616 

.663 

.922 

.6566 

.8735 

1.530 

3800 

4198 

6166 

1.62 

Porosity 

12 

34 

200 

293 

.637 

.686 

.918 

6566 

.8735 

1.530 

3800 

4624 

6826 

1.90 

16 

34 

25 

37 

.611 

.658 

.827 

.6004 

.8890 

1.4894 

3800 

3763 

5685 

1.50 

16 

34 

100 

147 

.617 

.664 

.864 

.6004 

.8890 

1.4894 

3800 

3985 

6030 

1.59 

16 

34 

200 

293 

.645 

.695 

.844 

.6004 

.8890 

1.4894 

3800 

4430 

6897 

1.82 

Ribbon 

12 

34 

25 

37 

.586 

.632 

.859 

.6558 

.8768 

1.5326 

3800 

3323 

5335 

1.40 

24%  Geometric 

12 

34 

100 

147 

.615 

.663 

837 

.6558 

.8768 

1.5326 

3800 

3714 

6163 

1.62 

Porosity 

12 

34 

200 

293 

.632 

.681 

.877 

.6558 

.8768 

1.5326 

3800 

4280 

6683 

1.76 

16 

34 

25 

37 

.603 

.650 

797 

.5570 

.8578 

1.4148 

3800 

3438 

5358 

1.41 

16 

34 

100 

147 

.626 

.674 

791 

.5570 

.8578 

1.4148 

3800 

3804 

5983 

1.57 

16 

34 

200 

293 

.648 

.698 

781 

.5570 

.8578 

1.4148 

3800 

4164 

6656 

1.75 

Cross  Chute 

34 

25 

37 

.710 

.543 

1.242 

.8867 

1.2776 

2.1643 

1928 

3768 

5798 

3.01 

W/L  =  .264 

34 

100 

147 

.707 

.540 

1.270 

.8867 

1.2776 

2.1643 

1928 

3810 

5712 

2.96 

34 

200 

293 

.716 

.547 

1.285 

.8867 

1.2776 

2.1643 

1928 

4212 

5925 

3.07 

47 

25 

37 

.759 

.580 

1.113 

.8494 

1.2512 

2.1006 

1928 

4052 

6868 

3.56 

47 

100 

147 

.729 

.557 

1.205 

.8494 

1.2512 

2.1006 

1928 

3973 

5958 

3.09 

47 

200 

293 

.775 

.592 

1.110 

.8494 

1.2512 

2.1006 

1928 

4292 

7303 

3.79 

REPRODUCED  FROM  REFERENCE  (6) 


TABLE  II  SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  24-GORE  AND  30-GORE  CONFIGURATIONS 


Parachute 

No.  of 

Suspension 

Velocity 

Scale  Factor,  K 

N 

Axes  Ratio 

Volume  in 

Type 

Gores 

Line  Length 

mph 

fps 

2a 

2a 

a 

b 

b' 

b  b' 

v 

V 

inches 

Do 

Df 

a 

a 

-  +  ~ 
a  a 

H 

c 

Flat  Circulars 

24 

34 

50 

73 

.677 

.679 

795 

.5758 

.8126 

1  3884 

4362 

4695 

30 

34 

17 

25 

.668 

.669 

.827 

.6214 

.7806 

1.4020 

4342 

4626 

10%  Extended* 

24 

34 

100 

147 

665 

.648 

.834 

.5949 

.8771 

1.4720 

4138 

4446 

Skirt 

30 

34 

17 

25 

.650 

.633 

.825 

.6255 

.7962 

1.4127 

4172 

4076 

Ring  Slot  16% 

24 

34 

25 

37 

.663 

.665 

.824 

.5800 

.9053 

1.4853 

3591 

3878 

Geometrically 

24 

34 

100 

147 

.680 

.682 

.819 

.5800 

.9053 

1.4853 

3591 

4079 

Porous 

24 

34 

200 

293 

.694 

.696 

.809 

.5800 

.9053 

1.4853 

3591 

4270 

30 

34 

25 

37 

.677 

.678 

.788 

.5800 

.9053 

1.4853 

3582 

3826 

30 

34 

100 

147 

.684 

.685 

.802 

.5800 

.9053 

1.4853 

3582 

4023 

30 

34 

200 

293 

.698 

.699 

.800 

.5800 

.9053 

1.4853 

3582 

4260 

Ribbon  24% 

24 

34 

25 

37 

.671 

.673 

.770 

.5980 

.8187 

1.4167 

3591 

3591 

Geometrically 

24 

34 

100 

147 

.676 

.678 

.813 

.5980 

.8187 

1.4167 

3591 

3927 

Porous 

24 

34 

200 

293 

.687 

.689 

.804 

.5980 

.8187 

1.4167 

3591 

4061 

30 

34 

25 

37 

.655 

.657 

.782 

.6021 

.8463 

1.4484 

3582 

3396 

30 

34 

100 

147 

.669 

.670 

784 

.6021 

.8463 

1.4484 

3582 

3622 

30 

34 

200 

293 

.677 

.679 

.823 

.6021 

.8463 

1.4484 

3582 

4002 

V 

o 


7273  1.67 

7027  1.62 


6930  1.67 

6265  1.50 


6031 

6510 

6924 

6404 

6588 

7012 


1.68 

1.81 

1.93 

1.79 

1.84 

1.96 


5968 

6097 

6389 

5666 

6022 

6256 


1.66 

1.70 

1.78 

1.58 

1.68 

1.75 


Since  this  parachute  was  “breathing"  during  the  test,  several  photographs  were  taken  at  each  speed.  The  data  were  reduced  from  the  photograph  which  most 
reasonably  appeared  to  represent  the  equilibrium  state. 


REPRODUCED  FROM  REFERENCE  (7) 
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X.  List  of  Symbols 

Ar.  -  Steady- state  projected  area  of  the 
Inflated  parachute,  ft2 

/\j  -  Instantaneous  canopy  mouth  area,  ft2 

A-,,  -  Steady-state  inflated  mouth  area, 

“  rtp 

a  -  Acceleration,  ft/sec^ 

-  Maximum  inflated  parachute  diameter 
of  gore  mainsoam,  ft 

b  -  Minor  axis  of  the  ellipse  bounded 

by  the  major  axis  (2a)  and  the  vent 
of  the  canopy,  ft 

b'  -  Minor  axis  of  the  ellipse  which 
includes  the  skirt  hem  of  the 
canopy,  ft 

C  -  Effective  porosity 

Cp  -  Parachute  coefficient  of  drag 

C  -  Parachute  pressure  coefficient, 

**  relates  internal  and  external 

pressure  (aP)  on  canopy  surface  to 
the  dynamic  pressure  of  the  free 
3  ire am 

U  -  Nominal  diameter  of  the  aerodynamic 
decelerator  =  ft 

F  -  Instantaneous  force,  lbs 

F  -  Siead3r-state  drag  force  that  would 

be  produced  by  a  fully  open  para¬ 
chute  at  velocity  Vs,  lbs 

F  -  Constructed  strength  of  the  para¬ 

chute,  lbs 

^rnax  ~  Maximum  opening-shock  force,  lbs 

o 

g  -  Gravitational  acceleration,  ft/sec 

k  -  Permeability  constant  of  canopy 

cloth 

m  -  Mass,  slugs 

M  -  Mass  ratio  -  ratio  of  the  mass  of 

the  retojrded  hardware  (including 
parachute)  to  a  mass  of  atmosphere 
contained  in  a  right  circular 
cylinder  of  length  (VstQ),  face  area 
(CpS0),  and  density  (o) 

M*  -  Mass  flow  ratio  -  zatio  of  atmo¬ 
sphere  flowing  through  a  unit  cloth 
area  to  the  atmosphere  flowing 
through  a  unit  inlet  area  at 
arbitrary  pressure 

n  -  Permeability  constant  of  canopy 
cloth 


P  -  Cloth  permeability  -  rate  of  air¬ 
flow  through  a  cloth  at  an  arbi¬ 
trary  differential  pressure, 
ft 3/ft 2/  sec 

p 

q  -  Dynamic  pressure,  lc/ff 

S  -  Instantaneous  inflated  canopy 

s  surface  area,  ft* 

9 

5  =A  -  Canopy  surface  area,  ft" 

o  so  9 

t  -  Instantaneous  time,  sec 

t  -  Reference  time  when  the  parachute 

has  reached  the  design  drag  area 
for  the  first  time,  sec 

tf  -  Canopy  inflation  time  when  the 
inflated  canopy  has  reached  its 
maximum  physical  size,  sec 

u  -  Air  velocity  through  cloth  in 
effective  porosity,  ft/sec 

v  -  Fictitious  theoretical  velocity 

used  in  effective  porosity,  ft/sec 

V  -  Instantaneous  system  velocity,  ft/sec 

V  ~  System  velocity  at  the  time  t  =  t  , 

°  ft/sec  ° 

Vs  -  System  velocity  at  the  end  of 

suspension  line  stretch,  ft/sec 

VQ  -  Volume  of  air  which  must  be  collec¬ 
ted  during  the  inflation  process, 

ft3 

W  -  Hardware  weight,  lb 

x^  -  Instantaneous  shock  factor 

X^  -  Shock  factor  at  the  time  t  =  t 
o  o 

p  -  Air  density,  slugs/ft^ 

r|  -  Patio  of  parachute  projected  mouth 

area  at  line  stretch  to  the  steady- 
state  projected  area 

g  -  Instantaneous  elongation 

emax  ”  ^ax^mum  elongation 

e  -  Initial  elongation  at  the  beginning 

of  the  elastic  phase  of  inflation 

S . F.  -  Parachute  safety  factor  =  Fc/Fmax 
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Appendix  B 

A  GUIDE  FOR  THE  USE  OF  APPENDIX  A 

/ 


At  first  reading.  Appendix  A  may  appear  to  be  a  complicated  sys 
tem  of  analysis  because  of  the  many  formulae  presented.  Actually, 
once  understood,  the  technique  is  straightforward  and  uncomplicated 
The  author  has  attempted  to  simplify  the  algebra  wherever  possible. 
This  appendix  presents,  in  semi-outline  form,  a  guide  to  the 
sequence  of  calculations  because  the  analysis  does  require  use  of 
formulae  from  the  text,  not  necessarily  in  the  order  in  which  they 
were  presented.  Also,  the  user  can  be  referred  to  graphs  of  per¬ 
formance  to  illustrate  effects. 

In  order  to  compute  tQ,  other  parameters  must  be  obtained  from 
various  sources. 

I.  Determine  System  Parameters 

1.  CdSq,  drag  area,  ft^  obtained  from  design  requirement. 

2.  V s ,  fps,  velocity  of  system  at  suspension  line  stretch. 

3.  p,  slugs/ft^,  air  density  at  deployment  altitude. 

4.  W,  lb,  system  weight  (including  weight  of  the  parachute) 
from  design  requirements. 

5.  V0 ,  ft^,  this  volume  of  air,  which  is  to  be  collected 
during  inflation,  is  calculated  from  the  steady-state  inflated 
shape  geometry  of  the  particular  parachute  type.  The  nomenclature 
is  described  in  Figure  24,p.A-14.  When  DQ  or  DF  is  known,  a  can  be 
calculated  from  data  in  Table  I  and  Table  II,  p.  A-15,  for  various 
parachute  types  and  number  of  gores.  Then  the  geometric  volume  VQ 
can  be  calculated  by  Equation  (31),  p.  a-14,  with  appropriate 
values  of  b/a  and  b'/a  from  the  tables. 

6.  AMq,  steady-state  canopy  mouth  area 
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where  N/a,  b/a,  and  b'/a  are  available  from  Tables  I  and  II  for  the 
particular  type  of  parachute  and  number  of  gores. 

7.  Aso ,  ft2,  canopy  surface  area  =  |do2 

8*  C  ,  pressure  coefficient,  see  Figure  18,  p.  a-12.  A 
constant  Cp  =  1.7  for  all  altitudes  seems  to  yield  acceptable 
results . 


9.  Constants  k  and  n  are  derived  from  measurements  of  the  air 
flow  through  the  cloth.  Only  k  is  needed  for  Equation  (14),  but  n 
is  also  required  for  Equation  (13) .  These  parameters  can  be  deter¬ 
mined  for  any  cloth  using  the  technique  described  beginning  on 
p.  A-12.  The  two-point  method  is  adequate  if  the  AP  across  the 
cloth  is  in  the  range  of  AP  for  actual  operation.  Check-points  of 
cloth  permeability  can  be  measured  and  compared  to  calculated 
values  to  verify  agreement.  If  the  data  are  to  be  extrapolated  to 
operational  AP ' s  greater  than  measured,  a  better  method  of  deter¬ 
mining  k  and  n  from  the  test  data  would  be  a  least  squares  fit 
through  many  data  points.  This  way  errors  due  to  reading  either 
of  the  two  points  are  minimized. 


II.  Step  1 


Calculate  the  reference  time  tQ  by  use  of  Equations  (13)  or 
(14),  p.  A-7.  If  the  deployment  altitude  is  50,000  feet  or  higher. 
Equation  (14)  is  preferred  due  to  its  simplicity.  For  altitudes 
from  sea  level  to  50,000  feet,  Equation  (13)  is  preferred.  Figure 
12,  p.  A- 8 ,  shows  the  effect  of  altitude  on  tQ  and  can  be  taken  as 
a  guide  for  the  user  to  decide  whether  to  use  Equation  (13)  or  (14). 
One  should  keep  in  mind  that  the  opening  shock  force  can  be  a 
strong  function  of  inflation  time,  so  be  as  realistic  as  possible. 

If  Equation  (13)  is  elected,  the  method  in  use  at  the  NSWC/WO  is 
to  program  Equation  (13)  to  compute  the  parachute  volume,  Vq,  for 
an  assumed  value  of  tQ.  Equation  (14) ,  because  of  its  simplicity, 
can  be  used  for  a  first  estimate  of  tQ  at  all  altitudes.  The  com¬ 
puted  canopy  volume  is  then  compared  to  the  canopy  volume  calculated 
from  the  geometry  of  the  parachute  as  per  Equation  (31),  p.  A-14. 

If  the  volume  computed  from  the  mass  flow  is  within  the  volume 
computed  from  the  geometry  within  plus  or  minus  a  specified  delta 
volume,  the  time  tQ  is  printed  out.  If  not  within  the  specified 
limits,  tQ  is  adjusted,  and  a  new  volume  calculated.  For  a  35-foot 
D0 ,  T-l 0  type  canopy,  I  use  plus  or  minus  10  cubic  feet  in  the 
volume  comparison.  The  limit  would  be  reduced  for  a  parachute  of 
smaller  D  . 

If  Y0  calculated  =  Y0  geometry  +  10,  then  print  answer. 


If  Y0  calculated  ^  Yn  geometry  +  10,  then  correct  tn  as  follows: 


Y  geometry 

Xo  *0  Y  calculated 
o 


(B-2) 
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The  new  value  of  tQ  is  substituted  in  the  "do  loop"  and  the 
volume  recomputed.  This  calculation  continues  until  the  required 
volume  is  within  the  specified  limits. 

III.  Calculate  tQ  corrected  for  initial  area.  The  tQ  of  Section 
II  assumes  that  the  parachute  inflated  from  a  zero  initial  area. 

If  this  is  a  reasonable  assumption  for  the  particular  system  under 
study,  then  the  mass  ratio  can  be  determined  from  Equation  (6) , 

p.  A-4.  For  h  =  0  if  the  value  of  M  <  0.19,  then  a  finite  state 
of  deployment  exists,  and  the  time  ratio  of  occurrence  and  the 
maximum  shock  factor  can  be  determined  from  Equations  (9)  and  (10), 
respectively,  on  p.  A-5.  If  p  ^  0 ,  then  the  limiting  mass  ratio 
for  finite  operations  will  rise  slightly  as  describee  in  Appendix 
C.  Figures  C-l  and  C-2  illustrate  the  effects  of  initial  area  on 
limiting  mass  ratios  and  shock  factors  respectively.  If  the  mass 
ratio  is  greater  th  a  n  the  limiting  mass  ratio  (ML)  ,  then  the 
maximum  shock  force  occurs  at  a  time  greater  than  tQ  and  the 
elasticity  of  the  materials  must  be  considered  (see  Section  VI) . 

If  n  ^0,  then  the  reference  time,  tQ,  will  be  reduced,  and 
the  mass  ratio  will  rise  due  to  partial  inflation  at  the  line 
stretch.  Figures  9  and  10,  p.  A-6,  illustrate  the  effects  of 
initial  area  on  the  velocities  and  shock  factor  during  the 
"unfolding"  inflation.  Equation  (15),  p.A-9,  can  be  used  to 
correct  tg  calculated  for  the  cases  where  n  =  Ai/Ac.  If  the  initial 
value  of  drag  area  is  known,  Equation  (16) ,  p.  A- 9,  can  be  used  to 
correct  tQ  and  rechecked  for  limiting  mass  ratios  versus  n  in 
Appendix  C. 

IV.  Opening  shock  calculations  in  the  elastic  phase  of  inflation. 

It  has  been  considered  that  from  time  t  =  0  to  t  =  tD  the  para¬ 
chute  has  been  inelastic.  At  the  time  t  =  tQ  the  applied  aero¬ 
dynamic  load  causes  the  materials  to  stretch  and  the  parachute 
canopy  increases  in  size.  The  increased  size  results  in  an 
increase  in  load,  which  causes  further  growth,  etc.  This  sequence 
of  events  continues  until  the  applied  forces  have  been  balanced  by 
the  strength  of  materials.  The  designer  must  insure  that  the 
constructed  strength  of  the  materials  is  sufficient  to  resist  the 
applied  loads  for  the  material  elongation  expected.  Use  of 
materials  of  low  elongation  should  result  in  lower  opening  shock 
forces  as  CDSmax  is  reduced. 

When  the  mass  ratio  of  the  system  is  greater  than  the  limiting 
mass  ratio,  the  elasticity  of  the  materials  and  material  strength 
determine  the  maximum  opening  shock  force.  The  maximum  elongation 
emax  an<^  the  ultimate  strength  of  the  materials  are  known  from 
tests  or  specifications.  The  technique  begins  on  p.  A-9. 

At  the  time  t  =  tQ,  calculate  the  following  quantities  for  the 
particular  values  of  M  and  p. 

a-  V0/Vs  from  Equation  (18),  p.  A-9. 
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b.  Xc  from  Equation  (22),  p.  A-10. 

c.  £q  from  Equation  (23),  p.  A-10. 

d.  Determine  CDSmax/CDS0  from  Figure  15,  p.  A-10. 

e'  aai1nUlate  thS  ■'■nf lati°n  time  ratio  tf/t0  from  Equation 

f.  Calculate  the  maximum  shock  factor  from  Equation  (19) 

p .  A- 9  .  '  ' 

g.  Calculate  the  opening  shock  force  Fmax.  =  *jFs  where 

FS  =  ^V|CDSo 

h.  Calculate  filling  time,  tf(sec) 

--(*) 

V.  In  order  to  simplify  the  required  effort,  the  work  sheets  of  Table 
B— 1  are  included  on  pages  B-5  through  B-9  to  aid  the  engineer  in 
systematizing  the  analysis.  The  work  sheets  should  be  reproduced  to 
provide  additional  copies. 
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Table  B-l .  Opening  Shock  Force 


lb. 

ft/sec2 

ft. 

slugs/ft2 

fps. 


ft. 


ft. 

ft.2 

ft.2 


ft.2 


ft.2 
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Table  B-l .  Opening  Shock  Force 

(cont 1 d) 

(3)  £max;  determine  maximum  elongations  from 
pull  test  data  of  joints,  seams,  lines,  etc.  Use 
minimum  £max  determined  from  tests. 

(4)  Cp;  pressure  coefficient 

h.  Steady  state  drag,  F  (lb),  =  ^pV^C^S 

°  o  2  w  Lj  O 


i.  Parachute  constructed  strength,  Fc  (lb);  deter¬ 
mined  from  data  on  efficiency  of  seams,  joints,  lines. 
Constructed  strength  is  the  minimum  load  required  to 
fail  a  member  times  the  number  of  members. 


3.  Force  calculations 

a.  Calculate  tQ  for  n  =  0 ;  eq.  14,  p.  A-7. 

r-  P9Vo  r 


t  = 


14W 


^vscDso 


CDSo 


2W  L  amo  '  Aso  k  /CpP 


-] 


Check  Figure  13,  p.  A-8,  for  advisability  of 
using  eq.  13,  p.  A-7. 

b.  If  T)  =  0,  proceed  with  steps  c  through  e. 

If  T]  41  0,  go  to  step  f. 

c.  Mass  ratio,  M;  eq.  6,  p.  A- 4 

2W 


M  = 


d.  If  M  <  4/21  for  r\  =  0,  then  finite  mass 
deployment  is  indicated. 

(1)  Time  ratio  at  X£  max;  eq.  9,  p.  A- 5 

1 


■(¥f 


(2)  Max  shock  factor,  x^;  eq. ,  10,  p.  A- 5 


_  J6_  /  21IVl\7 

C  49  y  4  y 


P 

Fc 


M 


t 

to  @Xj 


max 


xi 


max 


lb. 


lb. 


sec. 
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Table  B-l .  Opening  Shock  Force 

(Cont 1 d) 


(3)  Max  shock  force,  Fmax  (lb) 


F  =  X. 

max.  i 


e.  If  M  >  4/21;  then  intermediate  mass  or  infinite 
mass  deployment  is  indicated  and  the  elasticity  of 
materials  is  involved.  Calculate  the  trajectory  con¬ 
ditions  at  time  t  =  t  . 

(1)  Velocity  ratio  @  t  =  t  for  p  =  0 

V 


1  + 


7M 


(2)  Shock  factor  X  @  t  =  t  for  r\  =  0 


X  * 

o 


'  v 

_ o 

V0 


[’*m]  ,  - 

(3)  Initial  elongation,  eQ;  eq.  23,  p.  A-10 


X  F„ 

e  =  — °  s  e 
o  _  max 


Cr,S„ 


p.  A- 10 


(4)  Determine  —  max  from  Figure  15,  p.  A-10 

CDSo 

(5)  Calculate  inflation  time  ratio,  tf  ;  eq.  24, 


CDSmax 


(6)  Calculate  maximum  shock  factor,  x»  •  Pn 
19,  p.  A-9  -  —  1  max’  q‘ 


tf\6 


X: 


rvs  +  i 

17  X  ,1 

- 

V  7M 

o 

LV  7  J 

_ 

(7)  Calculate  maximum  shock  force,  Fmax  (lb). 


Fmax  FS 

max 
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% 


V 

V0 


CnS 


Damax 

cDso 


max 


max 
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Table  B-l .  Opening  Shock  Force 

(cont'd) 

(8)  Inflation  time,  sec  =  tf  =  tQ  (  —  ] 


f.  If  H  ^  0,  correct  tQ  for  initial  area  effects; 

cr  s  ■  V  ^ 


eq.  16,  p.  A-9 

t 

g.  Mass  Ratio,  M,  eq.  6,  p.  A-4 

2W 


I- 


**D  j 
CnS 

D  o 


M  = 


pgVst0CDSo 


t0  calculated 


h.  Calculate  limiting  mass  ratio,  MT 


M,  = 


L  3  (1  - 17) 


9  2_,  3  1 

u71  +m’+7 


If  M  <  Mt  ,  finite  mass  deployment  is  indicated  and  x. 

1  max 

can  be  determined  by  eq.  8,  p.  A- 5  by  assuming  values 
of  t/tQ  and  plotting  the  data  using  the  methods  of 
Appendix  C. 


i.  If  M  >  Ml,  then  intermediate  mass  or 
infinite  mass  deployment  is  indicated  and  the 
elasticity  of  materials  is  involved.  Calculate 
the  trajectory  conditions  at  time  t  =  t0. 


p.  A-9 


(1)  Velocity  ratio  @  t  =  tQ  for  n  +  0 ;  eq.  18, 


V 

_ o 

V0 


hh 


~T?)2  +  ail  zvl  +  „2 


(2)  Shock  factor  X  @  t  =  t  for  r|  ^  0;  eq. 
22,  p.  A- 10 


X  = 


1 


fv 


1  +  M 


(1  -T?)2  ai!  ~7?)  +  n2 

7  2 


2  \VC 


(3)  Initial  elongation,  £  ;  eq.  23,  p.  A-10 


XoFS 


CnS. 


(4)  Determine  —  max  from  Figure  15,  p.  A-10 


CDSo 
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Table  B-l.  Opening  Shock  Force  {Con 


(5)  Calculate  inflation  time  ratio,  _  ;  eq.  24, 

p.  A- 10  .  to 

,cA,x\i 

°DSo 


(6)  Calculate  maximum  shock  factor,  x^  max; 
eq.  19,  p.  A- 9  / t  \ 6 


(*y 


Xi 


(7)  Calculate  maximum  shock  force,  Fmax(lb) 


Fmax  =  Xj  FS 
max 


(8)  Calculate  inflation  time,  tf(sec) 
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